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Abstract 


The cold interstellar medium (ISM) plays a central role in the galaxy 
evolution process. It is the reservoir that fuels galaxy growth via star 
formation, the repository of material formed by these stars, and a sensi- 
tive tracer of internal and external processes that affect entire galaxies. 
Consequently, significant efforts have gone into systematic surveys of 
the cold ISM of the galaxies in the local Universe. This review dis- 
cusses the resulting network of scaling relations connecting the atomic 
and molecular gas masses of galaxies with their other global properties 
(stellar masses, morphologies, metallicities, star formation activity...), 
and their implications for our understanding of galaxy evolution. Key 
take-home messages are as follows: 

e From a gas perspective, there are three main factors that determine 
the star formation rate of a galaxy: the total mass of its cold ISM, how 
much of that gas is molecular, and the rate at which any molecular gas 
is converted into stars. All three of these factors vary systematically 
across the local galaxy population. 

e The shape and scatter of both the star formation main sequence and 
the mass-metallicity relation are deeply linked to the availability of 
atomic and molecular gas. 

e Future progress will come from expanding our exploration of scaling 
relations into new parameter space (in particular the regime of dwarf 
galaxies), better connecting the cold ISM of large samples of galaxies 
with the environment that feeds them (the circumgalactic medium in 
particular), and understanding the impact of these large scales on the 
efficiency of the star formation process on molecular cloud scales. 
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1. INTRODUCTION 
1.1. The nature and importance of the cold interstellar medium 


A detailed understanding of the interstellar medium (ISM) is crucial for most astrophysical 
pursuits. The ISM is shaped by pressure, turbulence and magnetic fields, consumed by 
star formation and supermassive black holes, heated and enriched by stellar winds and 
supernova explosions, and replenished by accretion from the Intergalactic Medium (IGM) 
or from galactic fountains (Shapiro & Fiela[1976]. Being at the crossroads of physical and 
chemical processes operating on such a range of scales, it is an incredibly rich source of 
information, but also a considerable challenge to model. 

The ISM spans in excess of six orders of magnitude in temperature and density, with 
pressure equilibrium dictating that it breaks down into three distinct phases 
[Cox][2005): the hot ionised medium (HIM), the warm medium, including 
both a neutral (WNM) and an ionized (WIM) component, and the cold neutral medium 
(CNM). In a galaxy like the Milky Way, the HIM and WIM dominate by volume, but due 
to their low densities (n ~ 107? — 10^! cm ?), most of the mass is in the cooler and denser 
neutral phases (Draine][2011). This cold ISM and its role in driving star formation and 
shaping galaxy evolution is the topic of this review. 

In a typical nearby star-forming late-type galaxy, the cold ISM accounts for ~ 10 — 20% 
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of the total baryonic mass (gas + stars). Of that cold ISM, ~ 75% by mass is atomic gas, 
mostly hydrogen (Hi), which is found in a warm (~ 5000 K), more diffuse (~ 0.6 cm?) 
phase and in a cold, denser phase (~ 100 K, ~ 30 cm ?), both contributing to the measured 
Hı emission at 21 cm (Draine[2011). The remaining ~ 25% by mass is in molecular form, 
with temperatures of ~ 10 — 100 K, and densities in the range of ~ 10? — 109 cm~?. Dust, 
amounting to ~ 1% of the gas mass, is mixed in with both the atomic and molecular 
components of the ISM, and therefore traces the total gas profile, as illustrated in Figure 


Star formation occurs once surface densities of ~ 10M pc ? are reached. Above this 


threshold, most of the ISM is molecular (Wong & Biitz[2002] [Bigiel et a1. [2008]. 'The picture 
for early-type galaxies is rather different, with the cold ISM only accounting for « 3% of 
the total baryonic mass. If present, any molecular gas tends to be centrally concentrated in 
a disc with a small scale length, but this gas can nonetheless have central surface densities 


similar to late-type galaxies 2013). The atomic gas of early-type galaxies 
exhibits a huge range of possible morphologies (Serra et al.12012). 
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Figure 1: Radial distribution of cold gas, dust and stars in “typical” late- and early-type 
galaxies of the nearby Universe. The spiral galaxy sketch is inspired by |Casasola et al. 
(2017) and|Bigiel & Blitz| (2012). The elliptical model is based on ATLAS?” and Herschel 


Reference Survey results (Serra et al.}2012 2013 2012). The dust 


curves have been scaled up by a typical gas-to-dust ratio for late- and early-type galaxies 


(De Looze et al [2020] Smith e sl .2012). 


Significant progress has been made in the past decade to turn observations of the mass 
and distribution of the cold ISM into a framework of scaling relations that connect the 
gas with the internal and external properties of galaxies, as well as their redshift evolution 
(cs. [Catinela et a 
2020). Observations of galaxies in the local Universe (loosely defined 
here as z « 0.05) have been of great importance in painting this picture; only there do 
current facilities allow us to measure both the atomic and molecular gas across large samples 
of galaxies representative of the entire population. This review is centered around these 
observations, and the insights they bring into key questions in galaxy evolution. 

While the central role of gas in galaxy evolution has long been recognized (e.g. 


1977 1980), these new gas scaling relations have triggered a resurgence of 
gas-centric galaxy evolution models (e.g.|Dekel et al./2009 2010 
2011) [Lilly et al.|/2013). These so-called “equilibrium” or “gas regulator” models rely on 
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Interstellar Medium 
(ISM): Interstellar 
matter within 
galaxies and filling 
the space between 
the stars. 


Intergalactic Medium 
(IGM): The diffuse 
hot plasma that fills 
the space between 
galaxies and traces 
the underlying dark 
matter cosmic web. 


Circumgalactic 
Medium (CGM): 
The hot ionised gas 
halo of galaxies 
located outside of 
the stellar 
component but 
within the virial 
radius. It forms the 
interface between 
the IGM and the 
ISM. 


Baryon cycle: The 
lifecycle by which 
gas is accreted onto 
galaxies from the 
IGM/CGM, feeds 
star formation, and 
is returned to the 
external 
environment due to 
nuclear and stellar 
activity. 


the balance between the various components of the baryon cycle: they equate the inflow 
rate of gas onto a galaxy (typically tied to the accretion rate of dark matter onto halos) to 
the rates at which this gas is locked up into stars via star formation and returned to the 
outside environment by outflows. Their success, especially given their simplicity, reinforces 
the view that gas is central to galaxy evolution. 

Numerical simulations, on the other hand, strive to implement all of the physics re- 
quired to produce the entire z ~ 0 galaxy population, not just those galaxies that are in a 
state of equilibrium. Large volume simulations however lack the resolution to capture the 
physics of gas cooling and star formation, instead relying on “subgrid” prescriptions. This 
includes for example the recipe for partitioning the cold ISM gas into H1 and Ha, as well 
as prescriptions for star formation and stellar feedback. Interestingly, modern simulation 
suites all reproduce the general stellar and star formation properties of z = 0 galaxies well, 
while relying on a wide range of such subgrid prescriptions. It is proving more challenging 
for these simulations to simultaneously reproduce observations of both atomic and molec- 
ular gas, but any discrepancies can be used to test and improve on the subgrid models 


for phenomena such as the Hi-to-H2 transition and stellar feedback (e.g. 2011 
[Crain et al.|2017 [Davé et aL.|2011 2020 [Diemer et al.|2018). 'Through these numerical 


efforts, global cold ISM measurements have the power to shed light on physics on much 
larger and smaller scales. 


1.2. Scope and objectives of this review 


This review is focused on how the total atomic and molecular gas masses of galaxies connect 
to their stellar mass, morphology, star formation activity, and other global properties, and 
how those dependencies fit in the broader galaxy evolution framework. Our main aims 
are to showcase how such global observations of the cold ISM can inform a range of key 
questions in the field of galaxy evolution (e.g. the nature and scatter of the galaxy main 
sequence, galaxy quenching, and the efficiency of the star formation process), by both being 
a direct agent in these processes, or by encoding information about phenomena that happen 
on much larger or smaller scales. 

Much about the link between the ISM and star formation has been learnt from carbon 
monoxide (CO) and HI observations at ~ kpc resolution of samples of up to ~ 50 — 100, 
mostly spiral galaxies (e.g. HERACLES; and EDGE; |Bolatto et al.|(2017) 
for molecular gas, THINGS; Walter et al.| (2008a) for Hr), and now even at resolution of 
tens of pc for CO studies (e.g. PHANGS; 2021). Results from such studies are 
touched upon as part of the discussions in Section [5] but will otherwise not be reviewed 


in detail; due to the unavoidable compromise between spatial resolution and sample size, 
they do not yet provide the full overview of the local galaxy population we are seeking 
here, but nonetheless are crucial in connecting global gas properties with relevant physical 
processes. Similarly, we do not review in detail the connection between the ISM and the 
larger scale gaseous environment, because an overall observational picture has yet to be 
fully established. This is an example of an area where there is significant scope for future 
discoveries with upcoming facilities, as discussed towards the end of this review. 
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2. OBSERVING THE COLD ISM 
2.1. Measuring the atomic ISM: the Hi 21-cm line 


Neutral atomic hydrogen (HI) is the most common element in the Universe, and the main 
constituent of the ISM of nearby disk galaxies. HI emits a spectral line at 21 cm, which was 
theoretically predicted by H. C. van de Hulst in 1944 and first detected in the Milky Way by 
(1951). 'The line arises from a magnetic dipole transition between the hy- 
perfine structure levels of the hydrogen atom in its ground state. The splitting of the ground 
state into two energy levels is caused by the magnetic interaction between the electron and 
proton spins (with the parallel spin state having higher energy), and the frequency of the 
electron spin-flip transition between the two is one of the most precisely measured physical 
quantities, vio = 1.420 405 751 786(30) GHz, corresponding to a wavelength of 21.106 cm. 
The transition probability is very small (Aio = 2.86888(7) x 10^ 1?s5^ 1), which translates 
into a radiative half-life for spontaneous de-excitation of 11/2 ~ 1/A1o œ 1.11 x 107 yr. In 
most astrophysical situations, the relative population of the hyperfine levels is determined 
by collisions. 

Because of its long wavelength, HI emission from external galaxies is spatially unresolved 
(except for the nearest systems) by single-dish telescopes, which have a spatial resolution 
(half-power beam width, HPBW) ranging from a few to tens of arcminutes. The construc- 
tion of radio interferometers such as the Westerbork Synthesis Radio Telescope (WSRT) 
in the Netherlands (completed in 1970) and the Very Large Array (VLA) in New Mexico 
(completed in 1980), opened up the possibility to map the spatial distribution and velocity 
field of the H1 emission within galaxies. However, the gain in spatial resolution of radio 
interferometers comes at a cost as their smaller collecting areas (i.e., lower sensitivity), 
compared to large single-dish telescopes, means that HI maps have become available almost 
exclusively for small samples of star-forming galaxies, generally within 100 Mpc from the 
Milky Way (c.g.,_van der Hulst et al.[2001| Walter et al.|2008b). Luckily, the next-generation 
Hi surveys on the Square Kilometre Array pathfinders (see Section [6) have already begun 
to increase these samples. 

Even without spatial information, single-dish global Hr-line profiles provide three impor- 
tant parameters: accurate redshifts (the spectral resolution of H1 observations of galaxies 
is typically a few km s ^! or better), total Hi-line fluxes and velocity widths (a measure of 
the Doppler broadening). If a galaxy is unresolved by the radio telescope beam, its total Hr 
mass can be simply obtained from its measured HI flux density integrated over the width 
of the line, assuming that the emission is optically thi] using the following expression: 


2.356 x 10? 


where Dz is the luminosity distance (in Mpc) to the galaxy at redshift z, as measured from 


the HI spectrum in the observed velocity frame, and the integral is the Hr-line flux density 


integrated over the line, with S, in units of Jy and v in km s !. 


1Much of the HI mass of galaxies is in a warm phase with spin temperatures of several thousand 
K, which is optically thin. This assumption does not hold in denser regions of the ISM (c.g., [Braun] 
Braun et al.[2009), but their filling factor is usually small. A correction for HI self-absorption, 
significant for galaxies close to edge-on view, has sometimes been applied to the Hr mass (see e.g. 


appendix B in|Haynes & Giovanelli]1984| and |Springob et al./2005). 
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The HI spectral line proved to be a powerful tool to investigate not only the atomic gas 
reservoirs of galaxies and their kinematics, but also to determine extragalactic distances 
and peculiar motions of galaxies (via the Tully-Fisher relation; and 
map the large-scale structure of the nearby Universe, especially in low-density regions of 
superclusters (where late-type galaxies are abundant and radio redshifts easy to measure). 
Generally observed in emission, but also in absorption against continuum radio sources, the 
21 cm line is not affected by interstellar dust and can therefore probe regions that are not 
accessible to optical telescopes. 


2.2. Measuring the molecular ISM: the CO molecule 


Most of the molecular ISM is composed of Hz. The molecule was first detected in the Galac- 


tic ISM via Lyman resonance absorption bands (Carruthers|}1970), confirming theoretical 


predictions that in regions of the ISM with visual extinction in excess of one magnitude, 
most of the hydrogen is in molecular form (Hollenbach et al.[1971). 

As for any other diatomic molecule, Hz has electronic, vibrational and rotational en- 
ergy levels. However, because H2 is homonuclear and therefore perfectly symmetric, the 
vibrational and rotational states can only emit radiation through their quadrupole moment. 
These lines have low transition probabilities, but even more problematically, high excitation 
energies requiring either gas temperatures in excess of 500K or a strong UV radiation field, 
conditions that are much more the exception than the norm when it comes to the molecular 
ISM. This leaves observers in the unfortunate situation that the species accounting for most 
of the mass of the molecular ISM is almost completely invisible to them. 

However, over 200 different molecules have now been detected in the ISM of the Milky 
Way and other galaxies (McGuire[2018). Since each transition of each of these molecules 
has different excitation requirements, observers are armed with invaluable tools to probe the 
temperature and density of the molecular ISM. Most of these molecules however have very 
low abundances compared to Hz, making them detectable in only the nearest or brightest 
of galaxies. For this reason, most of what we know about the molecular ISM of galaxies 
comes from observations of the second-most abundant molecule, carbon monoxide (CO). 
These observations are therefore the focus of this review. 

The lowest energy transition of the "CO molecule, at just hv/kp = 5.5 K and with a 
critical density of ~ 3000 cm ?, is the J = 1 — 0 rotational transition (hereafter referred 
to as the CO(1-0) emission lin’). These properties are conveniently very well matched to 
the typical conditions in the diffuse molecular phase of the ISM and in molecular clouds, 
which have temperatures of ~ 10 — 50 K and densities of ~ 10? — 10? cm? (Draine[2011). 
In addition, the wavelength of the CO(1-0) emission line, at 2.6mm, falls in a relatively 
transparent atmospheric window, making it accessible to ground based single-dish and 
interferometric millimetre telescopes. 

Most of what we know of the molecular ISM in galaxies of the nearby Universe therefore 
comes from observations of the CO(1-0) line. The observed CO(1-0) line flux is converted 
into a luminosity Loo (see e.g|Solomon et al.}1997), which is then itself turned into a total 


In this review, CO always refers to 12C190. Other species such as 1? C160 and 12C180 are also 
present in the ISM, but less widely observed as significantly fainter due to low abundances. The 
flip side is that while 1?C160 is generally optically thick, the isotopologues are optically thin and 
therefore more sensitive to the gas density. 
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molecular gas mass using the CO-to-H» conversion factor, aco: 
Mu, [Me] = acoLcoa—o) = aco [3.25 x 10 Sco vas Di (14-2) 7], 2. 


where Loo. oy is the line luminosity in units of K km s! pc”, Sco the integrated line flux 
in Jy km s E Vobs the observed frequency of the CO(1-0) line in GHz and Dz the luminosity 
distance in Mpc. For ISM conditions like the Milky Way’s, the conversion factor is generally 
taken to be aco= 3.2 Mo (K km s^! pc?) ! (this value is not inclusive of the correction 
to account for the contribution of helium and heavier elements; see Section [2.4}. 'There are 
comprehensive reviews of the use of CO line fluxes for measuring the molecular gas mass of 
galaxies in [Bolatto et al.] and [Tacconi et al] (2020), including a detailed account of 
why CO(1-0) *works" as a molecular gas tracer despite being generally optically thick, how 
aco is measured and how it may deviate from the Galactic value in certain environments. 
Indeed, a significant challenge with using CO(1-0) as our molecular gas tracer of choice 
is that aco is actually a conversion function, varying with the conditions of the ISM, in 
particular the dust abundance and the radiation field environment (e. g. [Israel]1997] [Accurso] 


et ai [017). 


2.3. Measuring the cold ISM: dust-based tracers 


It is also possible to infer the cold gas mass of galaxies from their dust contents, a method 
that has attracted significant interest in recent years with observatories like Herschel pro- 
viding a wealth of far-infrared (FIR) data. Assuming that dust is well mixed with both the 
cold atomic and molecular phases of the ISM, the total gas mass of a galaxy can be inferred 
from the dust mass, Maust, and a gas-to-dust ratio, dcp: 


Mgas = Mut + Mnp, pem Maust x dap. 3. 


The total dust masses required to apply this method are best inferred from observations 
that probe the peak of the FIR spectral energy distribution (100-200 jum) or even longer 
wavelengths on the Rayleigh-Jeans tail, where emission from cold dust dominates. They 


can be derived from the photometry using either empirical (e.g. |Scoville et al. a 


2014 2019) or physically-motivated models (e.g.|Draine & Li|2007). A 


detailed description of dust mass measurements techniques in nearby galaxies can E found 


in|Galliano et al.| Galliano et al.| (2018). The second ingredient in Equation []is the gas-to-dust mass ratio, 
dap. The value of dap is inked to metallicity (Z), and in many applications a simple scaling 


with metallicity of dgp«x Z ^ is assumed (Issa et al.|/1990} [Leroy et al.|/2011). Studies of 


the ócp(Z) relation however reveal important scatter across the local galaxy population, 
and a possible steepening of the relation at lower metallicity (Rémy-Ruyer et al.|2014] 
[Vis et al.|[2017). For nearby galaxies, the method summarised in Eq. |3|is best used to 
estimate total gas masses, Mgas, or to estimate My, when resolved maps of dust and HI 
are available. This allows for a pixel-by-pixel estimation of My,, after subtraction of the 
Myr contribution to the total gas mass. When resolved maps of the CO distribution are 
also available, the technique is also often used to estimate aco by combining Equations [2] 
and [3] ( e.g. [[szaei[1997] [Sandstrom ot al] Sandstrom et al.12013). 

An satan route is to use the extinction caused by the dust in the UV/optical 
(Bohlin et al.|1978). So far this has been explored mostly for molecular gas mass estimations, 
using the mE correlation between the Balmer decrement (or total V-band attenuation, 


Av) and Mu, (e.g. 2019| |Yesuf & Ho|2019). More accurate results can 
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be obtained from modeling that takes into account variations in the gas-to-dust and dust- 
to-metal ratios, possible deviations from the Case B recombination scenario, and various 
dust geometries (e.g. [Piotrowska et al.|[2020). When resolved 
spectroscopic observations are available, a correlation is found between Ay estimated from 
the Balmer decrement and both dust and Hə surface density maps, but with significant 
scatter and galaxy-to-galaxy variations 2020 
Investigations into the nature of any parameters that may be responsible for driving this 
scatter, including the contribution of HI, are necessary to enable accurate gas masses to be 
inferred. 

Care must however be taken in interpreting gas masses inferred from dust tracers, with 
specific calibrations performing better at reproducing either the atomic, molecular, or total 
cold gas masses of nearby galaxies (e.g. [Groves et al.]2015). For example, Janowiecki et al.] 
(2018) find that compared to other Herschel bands and more complex dust mass methods, 
the luminosity at 500um (L500) best predicts Mgas (0.15 dex scatter), while Liso most 
directly estimates My, (0.2 dex scatter). For relatively massive (M..Z; 10!? Mc) and metal- 
rich galaxies, FIR-predicted My, values are in good agreement with CO measurements (e.g. 
[Genzel et al.|[2015), with intrinsic scatter of « 0.2dex. However, the applicability of the 
method, whether using dust in emission or absorption, needs to be further investigated for 
quiescent galaxies, in the lower mass / metallicity regime, in the presence of an AGN, and 
as a function of physical scale. 


2.4. The contribution of helium and metals 


In the ISM, helium and heavier elements are mixed in with the hydrogen (be it atomic or 
molecular) and therefore contribute to the total mass. To account for this, it is standard in 
the molecular gas community to apply a 36% upward correction on the molecular hydrogen 
mass obtained from any of the above methods. For example, this is done by using a Galactic 
aco factor of 4.35 Mo[K km s^! pc?|^!, as opposed to the value of 3.2 Mo; [K km s^! pc?] ^! 
that would be needed to infer a “pure” molecular hydrogen mass (Bolatto et al.[2013]. This 
leads to the somewhat confusing situation where “ My,” is used to represent the total mass 
in the molecular ISM, when a symbol of “Myoigas” might be more appropriate. In contrast, 
the standard practice in the HI community is to use My; calculated using Equation 
accounting only for the hydrogen. The 36% correction is sometimes applied on a case-by- 
base basis, when a total mass for the atomic medium is sought. To ensure homogeneity and 
avoid confusion, we adopt the following convention across this review: 


e Quantities with the subscripts “Hr” refer to atomic hydrogen, without heavier element 
correction, as standard in the H1 21cm literature. 

e Quantities with the subscripts “H2” refer to molecular hydrogen only, without helium 
and heavier elements. This is a departure from common practice in the field, but 
makes for a more homogeneous presentation and easier comparison with the HI results. 
All quantities with the “H2” subscript can be scaled up by 36% for comparison with 
results in the literature when appropriate. 

e Quantities with the subscripts “mol” refer to the total molecular gas, including the 
contribution of helium and metals. 

e Quantities with the subscripts “gas” refer to the total cold ISM, therefore including 
Hi, H2, as well as helium and all the heavier elements. 
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2.5. Extragalactic cold ISM surveys 


Since the first extragalactic detections of HI in 1953 (Kerr & Hindman|1953 
1954) and CO two decades later (Rickard et al.|1975| [Solomon & de Zafra|1975 


1977), huge effort has gone into measuring the global amount, and more recently the 
detailed spatial distribution and kinematics, of these species in galaxies in the local Universe 
and beyond. Interestingly, the scope of these surveys has shifted from mapping the large- 
scale structure with the Hr-line (as 21 cm redshifts of star-forming galaxies can be obtained 
very quickly with single-dish telescopes), to investigations of the statistical properties of the 
cold gas content in large galaxy samples on one hand, and of the physics connecting cold 
gas and star formation in individual objects or small samples on the other. We focus here 
on Hi and Hə surveys that are particularly relevant to the study of global scaling relations, 
and refer the reader to[Giovanelli & Haynes] (2015), and references 
therein for more comprehensive historical overviews. 

Most of the global cold gas measurements available today come from blind Hr-selected 
surveys. The EM and most impactful of these, the HI Parkes All Sky Survey (HIPASS; 
Barnes et al.|2001) and the Arecibd?|Legacy Fast ALFA Survey (ALFALFA; Giovanelli et al. et Giovanelli et al. 
2005), emm over 5000 galaxies out to z Taal 0. "E ee ea ‘Meyer et al.]2004| Wong et al.12006) 
= ~31,500 sources out to z < 0.06 ( Caw p eeeetehs he surveys d 
been invaluable to enable statistical EE mE as 2 mass and velocity functions, as well 
as characterise the gas-rich galaxy population. The main drawbacks are their limited spatial 
resolutions (with a HPBW of ~15.5 arcmin for HIPASS and ~3.5 arcmin for ALFALFA at 
the rest frequency of the Hi line) and sensitivities, which lead to a bias towards the most 
gas-rich systems within their volumes. At the same time, targeted HI studies have focused 
on galaxy populations and/or environments that were not well probed by large-area blind 


HI surveys, such as early-type systems (e.g., Knapp et al./1985 2006 


and galaxies in different environments, from isolated to nearby groups and 
clusters (es, [Verdes-Montencgro et al 2001] [Kilborn et al-2009] [Chu et al [2009) 

The samples with molecular gas measurements are considerably smaller, due to the 
lack of blind CO-selected surveys at z ~ 0. Particularly remarkable was the FCRAO 
Extragalactic CO Survey (Young et al.|[1995], which measured CO(1-0) in 300 nearby 
galaxies and remained the largest comprehensive study of the molecular ISM in the local 
Universe for 25 years. Because of a strong correlation between infrared and CO luminosities 


(e.g. Sanders et al.||1991), a significant part of early extragalactic CO observations was 


um ie targeting luminous infrared galaxies (e.g. |Radford et al.||1991 


. Recognizing this bias toward “exceptional” galaxies (e.g. starbursts and interacting 
oA efforts were made to target a broader range of galaxies. Samples of galaxies 
with integrated molecular line observations grew to include “normal” non-interacting spiral 


galaxies [Braine et al.|[1993] [Sage|[1 , cluster spiral galaxies (Kenney & Young 
Bosell ot al [1997], early-type galaxies (Combes et al 007] [Krips et al PO10] [Young et al 
, galaxies with active nuclei (Helfer & Blitz 1993} [Sakamoto et al.[1999| 
2003 2011 


The past decade has seen a shift from targeting specific types of galaxies, to surveying 


, and isolated galaxies ( 


?'The 305-meter diameter Arecibo radio telescope tragically collapsed on December 1st, 2020, 
after nearly 60 years of operation that have led to numerous ground-breaking discoveries. Its 
contributions had a tremendous impact for HI astronomy, and in particular for the topic of this 
review. 


wunw.annualreviews.org * The Cold Interstellar Medium 


9 


Dark galaxy: A class 
of galaxies theorised 
to have large HI gas 
reservoirs yet no 
detectable stellar 
populations. 


Missing satellite 
problem: The 
under-abundance of 
observed low mass 
galaxies compared 
with the number of 
equivalent dark 
matter halos 
predicted by 
cosmological 
simulations. 
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This 
became possible also thanks to the availability of large multi-wavelength surveys such as 


SDSS (York et al.[2000), 2MASS (Skrutskie et al.12006), and GALEX (Martin et al.|2005), 


which measured stellar masses, star-formation rates and a plethora of other properties that 


samples representative of the entire galaxy population in both HI and CO lines. 


could be used to define representative samples for H1 and CO follow up. The largest single 
cold gas surveys in the nearby Universe are the extended GALEX Arecibo SDSS Survey 
(xGASS; 2018) and the extended CO Legacy Database 
for GASS (xCOLD GASS; [Saintonge et al||2011a[b] 2017), which obtained Hr and CO 
observations for 1179 and 532 galaxies, respectively, at 0.02 « z « 0.05. The objects were 
selected from SDSS and are representative of the full galaxy population with M.» 10? Mo. 
The strength of these surveys comes from their sample size, the homogeneity of all the 
data products, and the depth of the observations. As will be shown in Section [4] the latter 
is of particular importance as it makes it possible to infer the cold gas properties of even 
the most gas-poor sub-populations, via either spectral stacking or the setting of stringent 
upper limits. Other notable global cold gas surveys include the Hī census for the REsolved 
Spectroscopy Of a Local VolumE (RESOLVE; [Stark et_al.|[2016) survey, designed to be 
baryonic mass limited; the Herschel Reference Survey (HRS; Boselli et al.|2010), which 
includes 322 galaxies (225 of which have CO observations and 315 with HI data from the 
literature) that span the full range of environments, from low density regions to the core of 
the Virgo cluster; the ALLSMOG survey [Cicone et al.|2017), which 
observed CO(2-1) in 88 galaxies in the mass range 10°° «M. « 10'°Mo, increasing the 
coverage of parameter space towards lower mass galaxies. These surveys have been key 
in making cold gas a mainstream element in the multi-wavelength toolbox of extragalactic 
observers. 


3. GLOBAL STATISTICAL PROPERTIES 


Before delving into the details of the multi-parameter relations connecting the cold ISM 
masses of galaxies with their star formation activity and other global properties, we start 
by asking much simpler questions: how much molecular and atomic gas is there in the local 
Universe, and where is that gas predominantly located? 

In an attempt to resolve the “missing satellite problem", one of the science goals of 
the large-area blind H1 surveys described in Section [2.5] was to search for “dark galaxies". 
These objects would have significant HI gas reservoirs, yet be optically-dark. The existence 
of these galaxies would not only help to bridge the gap between the number of dark matter 
halos and the optically-bright observed galaxies, but also represent a challenge for star 
formation models; how can galaxies have so much gas yet be unable to form stars? It would 
also indicate that a complete census of the cold ISM contents of the local Universe would 
have to take into account all this gas not associated with optically-detected galaxies. Most 
of the candidate dark galaxies have turned out to be intergalactic gas clouds, the result of 
galactic interactions or environmental processes (Cannon et al.|2015); for example one of 
the first reported dark galaxies, VIRGOHI21 (Minchin et al.| 2005), was later found to be 
part of a large HI tidal feature stripped from NGC4254 by its fall into the Virgo cluster 
(Haynes et al.|2007). Given that massive clusters where such events can take place are 
rare, the contribution of such interstellar gas clouds to the total HI mass budget in the local 
Universe is negligible. 

While we cannot discount the existence of a population of very low mass gas-rich but 
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optically-faint objects below the detection limit of current blind Hr-selected surveys, these 
results nonetheless suggest that the vast majority of the cold atomic gas in the nearby 
Universe is located in galaxies with well-established stellar populations. We do not have 
similar blind molecular gas surveys at z ~ 0, but it is reasonable to expect that the above 
statement is even more true for molecular gas, given the minimum pressure and density 
requirements for molecules to form out of atomic gas. 


3.1. Cold gas mass functions 


We can now look at the distribution of this gas by drawing the atomic and molecular 
gas mass functions. Having established that the vast majority of the cold gas resides in 
the galaxy population traced by optical surveys, these mass functions can be built using 
either blind or targeted (mostly optically-selected) samples; once completeness effects have 
been taken into consideration they should produce similar results (Fletcher et al.|2021). 
For samples that are flux-limited, mass functions are most often built using the “1/Vmaz” 
method (Schmidt]1968). This method accounts for Malmquist bias by weighting each galaxy 
by the inverse of the volume in which galaxies of that specific brightness are detectable. 
Variations on this method exist to take into account and correct for large scale structure in 


the observed volume (Efstathiou et al.11988). 
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Figure 2: Left: The Hi mass function from the ALFALFA survey, and right, the molecular 
gas mass function from xCOLDGASS. Figures adapted with permission from [Jones et al.] 


(2018) and |Fletcher et al.| (2021). The parameters of the best-fitting Schechter functions 
(Eq. 


) are given in each case (where m. — log M*, the characteristic mass). 


Both the Ht and Hz mass functions (as well as the CO luminosity function) are well 


represented by a single (1976) function of the form: 


atl 


l 4. 


where M; would be either the HI or H2 mass (or their sum if interested in a total cold gas 
mass function). In this equation, M* is the characteristic mass at which the function tran- 
sitions from a power law with slope a to an exponential function, and $. is a normalisation 
factor. A Schechter function fit to the Hr and Hz mass functions is shown in Figure [2] 
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Most of the blind Hr-selected surveys described in Sec. [2.5] have been used to determine 
the Hr mass function. Differences between their results can be explained by the specific 
volume and depth of the observations. For example, shallow surveys tend to underestimate 
M* and the overall abundance of massive HI systems (being rare, they are easily missed 


if the volume observed is too small; [Martin et al.[2010). The footprint of the observations 


with respect to the large scale structure of the nearby Universe also has a direct impact on 
the Hr mass function; report a steeper low-mass slope a in the region 
around the Virgo cluster and a flatter slope in under-dense environments. This may be 
explained by a population of gas-rich galaxies in the filaments around clusters. Results for 
the shape of the HI mass function in specific galaxy groups can vary significantly, possibly 
highlighting the importance of the larger scale environment into which these groups are 


found (Stierwalt ot al-[2009] [Jones et a [2018). 


Measurements of the H2 mass function are far fewer between due to the challenges in 
assembling the necessary samples. built a CO luminosity function and H2 
mass function from the FCRAO Extragalactic CO Survey (Young et al.}1995), for both a far 
infrared flux-limited sub-sample, and an optical B-band-selected sub-sample. To avoid the 
possible biases associated with the IR- or optical-flux selection, robust Hə mass functions 
have recently been computed using the large HRS and xCOLD GASS mass-selected samples 
(Andreani et al.|2020| [Fletcher et al.[2021). In addition to the sample selection, the choice of 
the CO-to-H» conversion function also has a significant impact on the H2 mass function and 
the best-fit Schechter function parameters 
poz). 


3.2. How much HI and Hg is there in the nearby Universe? 


By integrating the mass functions, we obtain a measurement of Qur and Qy,, the cosmic 
mass densities of H1 and Hz in the nearby Universe. Amongst many other measurements 
of Qur (e.g. [Zwaan et al.|2005), the ALFALFA survey estimates that Qui= (3.9 + 0.1 + 
0.6) x 1074 (Jones et al. 2018), while for molecular gas |Fletcher et al.| (2021) estimate 
Qmot = (7.62 + 0.47) x 107% from the xCOLD GASS sample, and 


report a value of (1.1 + 0.4) x 1074 for the HRS sample. Once accounting for the presence 


of Helium and heavier elements by multiplying Quı by a factor of 1.36 (already included 
in Qmo1), the total density parameter for the ISM in the nearby Universe is found to be 
QOrsm = 0.00061. Accounting for dust would increase this number by ~1% given standard 


dust-to-gas ratios measured in the local Universe (e.g. Rémy-Ruyer et al||2014). This 


corresponds to ~35% of the total mass density of stars, or ~1.3% of all the baryond?] 


3.3. How is the cold interstellar medium distributed across the nearby galaxy 
population? 

The good agreement between the H1 and Hz mass functions and a single Schechter function is 

in contrast with the stellar mass function in the nearby Universe, which is best characterised 


by a double Schechter function (Baldry et al.|2012). The double Schechter function, being 


the sum of two functions with the same M" but different a and ¢. values (as in Eq. 14): 
better reproduces the total stellar mass function by separately accounting for the star- 


4These values assume that Q. = 0.0017 (Baldry et al./2012) and Q = 0.0456 (Planck Collabo- 


2020), with all values calculated for Ho = 70 km s-! Mpc-!. 
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forming galaxies (abundant at low masses, with a high value for a and lower à.) and the 
quiescent population (dominating at higher masses, with low a and higher ¢.). If confirmed 
by more robust statistics at the low mass end, the fact that HI and Hz mass functions are 
well fitted by a single Schechter function suggests that, unlike stellar mass, most of the 
cold ISM mass in the nearby Universe is contained in a single galaxy population, namely 
star-forming galaxies as will be shown in Sec. 

The values of Qur and Qu, imply that, overall, the molecular-to-atomic gas abundance 
ratio in the local Universe is 14.4%. This is in stark contrast with the situation at z ~ 1.5, 
the point in cosmic history when molecular gas was at its most abundant with Qu, /Qur~ 
75% (Walter ct al[2020). 

There is also a stellar-mass dependence to the overall distribution of Ht and Hə in the 
local Universe, with low mass galaxies (M..« 10? Mo) contributing between 30% and 40% 
of the overall Qur, but only 11% of Qg,. Conversely, high mass galaxies (M.> 10/9 Mc) 
account for ~ 65% of the molecular gas, and ~40% of all the Hī in the local Universe 
[Fletcher etal 2021). Such 
detailed measurements are proving very valuable to test and constrain cosmological hydro- 
dynamic simulations, which while overall producing qualitatevely accurate HI and H5 mass 
functions and scaling relations, sometimes over- or under-predict the relative abundance of 
the two at either low or high stellar masses (see for a detailed discussion 
of these matters for the SIMBA, Eagle and IllustrisTNG simulations). 


4. COLD GAS SCALING RELATIONS 


We now turn our attention to the relations between global cold gas content and other galaxy 
properties, such as stellar mass, colour or star-formation rate — these are generally referred 
to as cold gas scaling relations, and usually expressed in terms of gas mass (Mar, Mu, or 
Mgas) or gas-to-stellar mass ratio (far= Mmur/M., fu; Mm,/ M. or fgas = Mgas/ Mx, 
often improperly called “gas fractions"). Because galaxy properties scale with mass, the 
latter formulation is usually preferred. 

The usefulness of gas scaling relations was recognised early on, as a means of quantifying 
the typical H1 and H» gas content of galaxies and deviations from it 
Young et al.|1995). In these studies, gas mass was typically normalised by luminosity 
in B-band or dynamical mass, and plotted as a function of luminosity, morphological type, 
colour, or other optical quantities. 

The study of gas scaling relations has undergone a renaissance in the past decade, thanks 
to the availability of large, representative samples of galaxies with HI and H2 measurements, 
along with a wealth of accompanying multi-wavelength data. Modern studies parameterise 
gas content in terms of stellar mass, stellar surface density, star-formation rate and other 
quantities that have an immediate physical interpretation and usefulness for theoretical 
models, facilitating the comparison with numerical simulations and semi-analytic models 
of galaxy formation and evolution. We review in Section [4-1] the main approaches used to 
compute scaling relations in the presence of upper (or lower) limits and missing data. 


4.1. Methods 


There are several aspects to keep in mind when it comes to quantifying gas scaling relations 
and comparing the results between different works. We discuss these below and refer to 
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Figure [B] which shows HI gas fraction as a function of stellar mass computed with different 
samples and/or methods to illustrate our main points. 


1.5 
€ ~ log<fHI>, A100 
@ stacking, A100 
1.0 @  «log(fHI)», XG 
E log<fHi>, xG 
@  median[fHI], xG 
0.5 * survival, xG 
0.0 
S 
3 -0.5 
—1.0 
-1.5 
-2.0 
9.0 9.5 10.0 10.5 11.0 11.5 
log(Mx)[Mo] 


Figure 3: Relation between HI gas fraction and stellar mass for HI detections (blue circles) 
from the ALFALFA survey (or “A100”, [Haynes et al.| (2018), as referred to in the legend) 
and xGASS galaxies (gray, with downward arrows indicating non-detections). The scaling 
relation between M.-f#r is derived in bins of M. using different methods: average fxr of 
ALFALFA detections (blue), spectral stacking of ALFALFA data (red; mostly individual 
non-detections), for xGASS the difference between the logarithm of linear averages (orange) 
and average of logarithms (yellow), the median value of logfg; (green; including both 
xGASS detections and non-detections), and the logarithmic averages from survival analysis 


(purple). 


Sensitivity of observations. Naturally, the sensitivity of the surveys used to compute 
scaling relations matters. For instance, unless one is interested in quantifying the average 
gas content of star-forming galaxies only, obtaining scaling relations from relatively shallow 
observations (such as wide-area blind HI surveys) using only detections will return results 
that are biased toward high gas content, compared to samples that span the full range of 
galaxy properties, from star-forming to passive systems. This is shown by the blue line 
inF ig.[3] obtained by computing (logarithmic) averages of ALFALFA detections (blue dots). 


Averaging procedure. There are several common estimators of the typical gas content, 
such as median and linear or logarithmic mean, which generally lead to different results, 
depending on the underlying distribution of the data. Hence, linear averages are usually 
offset from medians or logarithmic averages, reflecting the fact that galaxy properties 
typically follow log-normal distributions, and sometimes not even log-normal (see e.g. the 
bimodality of SFR in local galaxies). In Fig.|3| we use xGASS detections and upper limits 
(gray dots and arrows, respectively) to show the difference between logarithms of linear 
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averages (orange line), logarithmic averages (yellow line) and medians (green line). 


Treatment of non-detections and upper limits. Representative cold gas samples 
are typically selected from optical properties and followed up with radio and millimetre 
telescopes to measure the H1 and CO content, but not all galaxies are usually detected. 
Nonetheless, upper limits carry useful information that needs to be taken into account in 
the scaling relations, and there are three main approaches to do so. 

'The simplest option is to compute medians, rather than averages, to minimise the 
impact of the non-detections. Indeed, as long as the detection fractions are above 5096 
and the upper limits do not significantly overlap with the detections (as is the case for the 
xGASS and xCOLD GASS surveys, which were designed to reach gas fraction limits of a 
few percent to produce stringent upper limits), medians are unaffected by the exact value 
of the gas mass of non-detections. 

Another technique to deal with upper limits is spectral stacking, which has become 
a common tool to constrain the statistical properties of a population of galaxies that lack 
individual detections in a survey. This turned out to be very powerful when applied to 
wide-area, blind HI surveys, provided that independent spectroscopic redshifts were also 
available (e.g [Meyer ef al|B016). Brief 
the Hr spectra of N galaxies are aligned in redshift and co-added (or stacked), regardless 
of whether these are individually detected or not, resulting in a total spectrum that yields 
an estimate of the total and average (when divided by N) HI content of the galaxies that 
were co-added. The red line in Fig. [3] shows the result obtained by stacking the ALFALFA 
HI spectra from [Brown et al.[2015] However, this technique has two important limitations. 
First, spectral stacking is an intrinsically linear operation, hence scaling relations from 
stacking cannot be blindly compared to ones expressed in terms of averages of logarithmic 
quantities (e.g., compare the red and orange lines in the figure). This means that spectral 
stacking is most useful when used to compare relative offsets between subsets of galaxies 
(e.g., divided by specific SFR or halo mass; rather than absolute 
values. Second, a stacked spectrum carries no information on the scatter of the underlying 
population whose spectra were co-added; while a handle on the variance of the result is 
usually obtained via the so-called jackknife statistics or by binning the sample by multiple 
parameters simultaneously (if the size of the sample allows it), this remains challenging. 

A third strategy to take non-detections into account is to adopt survival analysis 
methods, which allow to compute linear regressions or means, medians and other statistical 
estimators for samples that include censored data such as upper limits, which are referred to 
as "left-censored data" in survival analysis. We refer the reader to (1985) 
for an introduction of this branch of statistics to astronomers, MI EE I 
discussion of existing approaches to compute likelihood functions in the presence of censored 
and missing data. These techniques have been applied to both Hr and Hə datasets (e.g., 
[Calette et al.[2018} 2021] [Feldmann[2020). Following|Feldmann] (2019), we show the averages 
for xGASS using survival analysis as a purple line in Fig.|3| The comparison with median 
values (green line) shows an offset of ~0.2 dex between the two methods, especially at the 
low-mass end where upper limits are fewer. This is just a consequence of the fact that the 
underlying gas fraction distribution is not log-normal, hence medians and averages do not 
give the same answer. 

Which technique is most appropriate depends on the sample used and on the specific 
question that one is trying to answer. Medians are meaningful only when upper limits are 
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well separated from the detections and do not dominate the statistics in the bin. When 
upper limits and detections are mixed, survival analysis should be the preferred approach 
(see e.g. [Stark et al.|/2021), considering that spectral stacking is intrinsically linear and 
thus should be used only to compare relative offsets rather absolute values. When bins 
are dominated by upper limits, determining scaling relations becomes more challenging, 
although survival analysis can recover means for detection fractions as small as 2096 (but 
with a large uncertainty; [Calette et al.]2018). In this review, we compare medians to 
averages from survival analysis to demonstrate the impact of the non-detections on the 
observed trends. 


4.2. Gas fraction scaling relations 


The Hı and H2 masses of galaxies correlate with several of their other global properties, 


such as luminosity, mass, size, morphology, and star formation rate (e.g. Kenney & Young 
1989 1993| Roberts & Haynes||1994). To illustrate our discussion of these results, 


Figures |4| and |5| show the main cold gas scaling relations from the xGASS and xCOLD 


GASS datasets as initially presented in |Catinella et al.|(2018) and (2017) 
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Figure 4: Gas mass fraction scaling relations between the atomic gas mass fraction 
(far=Mutr/M., top row), the molecular gas mass fraction (fH,=Mx,/M., middle row) 
and the molecular ratio (My, /Mg;, bottom row) as a function of stellar mass (M,, first 
column), and stellar mass surface density (u+, second column). The gray dots and limits 
are the xGASS/xCOLD GASS galaxies, and the binned values are calculated as medians 
(green hexagons) and mean values using survival analysis (purple crosses) as described in 


Sec. 


5While including other samples in these figures would extend the dynamic range in one or more 
dimensions (e.g., probing lower stellar masses, more gas-rich systems, interacting galaxies, early- 
type morphologies and so on), this would complicate the simple selection function of xGASS and 


xCOLD GASS and therefore the interpretation of any trends. See however (2018) 
and|Ginolfi et al.|(2020) for more extensive data compilations. 
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4.2.1. Correlations with structural properties. Figure shows the HI gas fraction 
(fur: Mur/M.; top row), Ha gas fraction (fu, Mr, /M.; middle row) and the molecular- 
to-atomic gas mass ratio (bottom row). These quantities are correlated with two basic 
structural properties: the stellar mass (first column) and the stellar mass surface density 
(p M.(2529,,) !). This quantity is commonly used as a quantitative measurement of the 
morphology of galaxies, with bulges growing in prominence as u+ increases and becoming 
dominant at log(u./ Mcokpc ?)» 8.7. 

Both far and fp, are anti-correlated with M+. Part of this trend is due to the increasing 
fraction of gas-poor, quiescent and mostly early-type galaxies as mass increases, although 
the trend is also found when considering only star-forming galaxies (Cicone et al.|[2017). 
There is also an anti-correlation between fg; and fg, and us, with galaxies becoming 
progressively more gas-poor the more bulge-dominated they are. Similar results are found 


when visual morphologies are used (e.g.,|Boselli et al./2014b). 


Of note is the flattening of both the M.-fz, (also noted by 2015) and 


[ix-fH, relations (Fig. [ap&a) below the characteristic thresholds of log(M..)e- 10.3 and 
log(j4.)e 8.5 where the galaxy population transitions from being dominated by quiescent, 
mostly early-type objects, to star-forming galaxies (Kauffmann et al.|[2003). The lack of 
strong dependence of molecular gas contents on M. and yx in late-type systems (which 
dominate in Fig. |4|at low M. and u+) has long been reported (e.g. 
2014b). 'This break is however not seen in the equivalent 
fur scaling relations (Fig. |4pé&b). indeed showed that the My,/Mur ratio 
increases amongst spiral galaxies as the bulge component grows in importance (as in Fig. 
[4f). They however pointed out that this result was possibly of little physical significance 
since the excess H1 in low-mass, disc-dominated galaxies is located at large radii, away from 


the regions where most of the molecular gas is found — we come back to this point in Section 
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Figure 5: Gas mass fraction scaling relations as in Fig. [4] but for NUV-r color and specific 


star formation rate (SSFR). The dashed lines in panels a and d correspond to Equations [6] 
and [5] respectively. 
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4.2.2. Correlations with star formation activity. In Figure[5] the gas quantities are shown as 
a function of two different estimates of the present-to-past averaged star formation activity, 
the NUV-r color and the specific star formation rate, SSFR. Both quantities compare a 
measurement of the current star formation activity (NUV flux in one case, the SFR in the 
other), with a measurement of the older stellar population (r-band flux and total stellar 
mass, respectively). The difference is in the modeling assumptions that have gone into the 
computation of SFR and M,, including a correction for dust extinction. 

The tightest of all the relations in Fig. Blis the one between fg, and SSFR. This well- 
known correlation between molecular gas and star formation activity is often also portrayed 


as the relation between CO and total infrared luminosities (e.g. Sanders & Mirabelļ||1985 


Solomon et al.]1997) or as the relation between the surface densities of molecular gas mass 
and star formation rate (ie. the Kennicutt-Schmidt relation, e.g. 
[et a1. [2008). In contrast, fg; correlates more strongly with NUV-r than with SSFR. One 
might intuitively expect that HI, being a long-term gas reservoir, would correlate more 
with a longer-term measure of SFR (i.e. traced by NUV) than a shorter-term one (e.g., 
traced by Ha or FUV). However, the SSFRs in Fig. Blare based on NUV-r (and infrared) 
measurements, so this is not the case. Instead, this can be explained by the correlation 
between HI and un-obscured star formation, which is evident in the outer parts of the star- 
forming disks, as illustrated convincingly by the case of M83 (Bigiel et_al.][2010}. More 
interesting than the correlation between gas and star formation itself however is the scatter 
of these relations, which reveals differences in global star formation efficiency across the 
galaxy population. This is discussed in Section 


4.3. Depletion time scaling relations 


Once it has been established how much atomic and molecular gas is typically found in 
galaxies of a given type, we can investigate how the presence of that gas impacts on their 
star formation activity. This is most commonly studied through the global star formation 
efficiency, SFE, defined as the ratio between the SFR and the gas mass (either atomic, 
molecular, or total), or conversely the depletion time, taep = (SFE) !. The depletion 
time is an indication of the timescale over which star formation could be maintained at 
the current rate, given the available gas supply, and assuming no inflows of fresh fuel or 
recycling. In the nearby Universe, the typical depletion time for molecular gas in star- 
forming spiral galaxies is ~ 1 Gyr, and this value decreases slightly towards higher redshifts 
[Saintonge et al.|2013). These short depletion times (« 1 Gyr) at 
look-back times of > 10 Gyr are used as evidence for the steady accretion of gas onto 
galaxies in the intervening time, or at least for transport of atomic gas from the outer discs 
imwards; without such re-supply, there would be no gas and therefore no star formation in 
the Universe today, clearly at odds with observations. Another limitation of this simple 
interpretation of tgep is that not all the gas present can participate in the star formation 
process; this is particularly true of HI, most of which does not have the density required 
for star formation and will remain as such in the absence of secular or external processes 
to change the status quo. 

Figure [6] shows how the atomic and molecular gas depletion times vary as a function 
of the same four parameters used in Figures |4| and Once upper limits are taken into 
account, taep(HI) shows no significant dependence on any of the quantities as judged by 
the Kendall 7 coefficients, but rather is characterised by a large scatter of 0.5 dex. The Ht 
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Figure 6: Atomic and molecular gas mass depletion time scaling relations, with parame- 
ters, data, and symbols as defined in Fig. In each panel, the Kendall rank correlation 
coefficient (7) is given. 


distribution typically extends well beyond the optical disc, meaning that significant (but 
variable) amounts of the total HI reservoir are not directly involved in the current star 
formation process. 

The molecular gas depletion time, t4.5 (H2), shows a different behaviour; it correlates 
with global galaxy properties (especially those with a direct link to star formation, like 
SSFR and NUV-r), and shows far less scatter. This is unsurprising given the much closer 
link between molecular gas and star formation, but it is a departure from the idea of a 
universal relation between molecular gas and star formation (e.g. [Bigiel et al. [2011], and a 
result that only emerged when parameter space was expanded sufficiently by surveys such 
as x COLD GASS for these trends to appear. The results in Figure |6]suggest instead that a 
given amount of molecular gas does not always produce the same levels of SFR; rather, the 
gas-to-star conversion efficiency appears to vary systematically across the galaxy population. 
Once combined with the variations of far and fx,, the implication of Figure [6] is that the 
specific level of star formation activity in galaxies is determined by both the amount of gas 
available, and the efficiency of the conversion of that gas into stars (e.g. [Saintonge et al. 
poro) 

The dependence of tae; (H2) on global galaxy properties is at first puzzling; how can the 
details of star formation, a process occurring on parsec scales, be “controlled” by the global 
properties of galaxies on scales of tens of kpc? Or could this apparent correlation be caused 
by another common factor? These questions are best answered by looking at the Kennicutt- 
Schmidt (KS) relation. On global scales, the taep(H2) variations shown in Fig. [6] manifest 
themselves as systematic deviations of particular galaxies. For example, bulge-dominated 
galaxies (e.g. those with high values of jz. in Fig. |6) tend to fall systematically below the 
KS relation traced by spiral galaxies (Saintonge et al.2011b} Davis et al.[2014). At ~ kpc 
resolution, the variations in taep(H2) are found to be larger from galaxy to galaxy than 
within individual galaxies; in other words, individual galaxies follow their own KS relation, 
but the exact slope and normalisation of this relation varies systematically based on global 
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galaxy properties (Ellison et al.|2021b). Variations are even seen at the Giant Molecular 


Cloud (GMC) level. In nearby galaxies where molecular gas has been mapped on scales of 
tens of pc, GMCs are found to be systematically bigger/brighter in some galaxies 
[et al.[2013a}. What’s more, these galaxy-to-galaxy variations in cloud properties are linked 
to global galaxy properties, such as mass and morphology, via changes in internal pressure 
and structure and therefore systematic variations 
in the amount of dense star-forming gas compared to the total (diffuse) molecular gas mass 
traced by CO(1-0). This is discussed further in Section [5.5] 


4.4. Application of scaling relations: estimating cold gas masses 


Gas fraction scaling relations, like those shown in Figures |4| and |5| give us expectation 
values for the gas contents of “normal” galaxies, which is of crucial importance to study 
the impact of additional factors such as presence of an AGN, interactions, or environment. 
For example, the tight correlation between Mg; and optical isophotal diameter at fixed 
morphological type was of particular importance early on as 
it led to the definition of the Hr deficiency parameter to quantify environmental effects 
and show how clusters impact on the ISM of infalling galaxies 
Cortese et al. [2021]. 


Stellar mass by itself is a relatively poor predictor of gas content, with a spread of over 
2 orders of magnitude in fy, and fxr at fixed M. (Fig. [4pé&c). This is especially true for 
high mass galaxies (M.> 10'°Mo), and unless additional information can be used to first 
break the degeneracy between star-forming and quiescent galaxies at fixed mass 


[et al. [2015]. There are however a number of strategies at our disposal: 


4.4.1. Control samples. The first approach is to use a large representative galaxy sample 
with direct Hi, CO (or dust) measurements, and to draw control samples from it. The Hr 
and H3 masses of a galaxy of interest can be estimated to be the mean of the gas masses of 
the galaxies most similar to it in the reference sample. The choice of which parameters to 
control against depends on the specific science goals and the properties of the sample, but 
the general gas scaling relations can guide those choices. Given all the correlations presented 
in Figures[4]and[5} it is recommended to select control galaxies by matching with the target 
object on two parameters, one that controls for mass (either M+, or luminosity in a NIR 
band) or alternatively morphology (u+ for example), as well as a parameter that traces star 
formation activity (SFR itself, or luminosity in the UV or FIR). The combination of M, 
and SFR is most often used to select control samples to establish how specific phenomena 


impact the gas contents of galaxies (e.g. 2018} |Ellison et al.|2019 
2021). 


4.4.2. Empirical relations. Alternatively, empirical scaling relations can be used to produce 
expectation values of My; and My,. These methods are best used in the same spirit as 
the control sample technique described above. If empirical estimates of Mg; and My, 
are nonetheless used in lieu of direct measurements (see Section [2.3), caution should be 
taken (1) to take into account the intrinsic scatter of the relation used, (2) to stay within 
the parameter space where the relations have been calibrated, and (3) to be mindful of 
underlying assumptions in the calibration, which may limit the applications. 
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4.4.2.1. Inferring My, from the SFR. Of all the scaling relations shown in Figures 
those directly relating My, and SFR are the tightest. This strong correlation has 
long been exploited to infer My, values, a technique sometimes referred to as “inverting the 
Kennicutt-Schmidt relation”. Based on the xCOLD GASS data (Fig. |5H) the recommended 
relation to predict My, values for galaxies in the local Universe is: 


log( Mg, /M.) = (0.75 + 0.01) log(SSFR) + (6.24 + 0.12). 5. 


A factor of log,)(1.36) should be added to the right hand side to obtain the total molecular 
mass fraction, log(Mmoi/M.), including the contribution of He and metals. This relation 
is best applied for galaxies with log(SSFR) » —11.5 where it is robustly constrained by 
observations. For CO detections, the scatter around this relation is 0.24 dex, of which 0.12 
dex is intrinsic, and 0.20 comes from measurement uncertainties. By combining local and 


high redshift observations, (2020) propose a calibration of this relation that 


is applicable at all redshifts up to z ~ 3. Others have fitted the 3D parameter space formed 
by Mp, M. and SFR to provide redshift-independent gas predictors (Santini et al. [2014]. 

For local Universe galaxies, the good availability of mid-infrared photometry from the 
WISE all-sky survey at wavelengths of 3.4, 4.6, 12, and 22 um can also 
be exploited to estimate My, values. The two shorter wavelength bands trace the estab- 
lished stellar population, while the 12, and 22 wm bands probe warm dust and Polycyclic 
Aromatic Hydrocarbon (PAH) emission, therefore being sensitive to the ongoing rate of star 


formation. (2017) have indeed shown that the flux ratio between 12 and 4.6um 


correlates well with fg,. This is explained by the fact that the 12-to-4.6jm flux ratio is a 
proxy for SSFR, which itself has been shown to correlate tightly with fr (see Section |4). 
A tight correlation (rs ~ 0.9) is also found directly between My, and Lizum 
2015). This is expected for star-forming galaxies where ~ 80% of the 121m emission is due 


to stellar populations younger than 600 Myr (Donoso et al./2012), but caution must be ap- 


plied in quiescent galaxies where older stellar populations contribute significantly, breaking 
the relation between MIR emission and molecular gas (Davis et al.[2014). 
expand on this by showing that adding M. and g — r to Lizum, yields predicted Lco(1.0) 
values with 0.16 dex of intrinsic scatter compared to direct measurements, an improvement 
from the 0.21 dex of scatter obtained by using Lizum only. 

When estimating My, values from the SFR, the implicit assumption of a constant 
molecular gas depletion timescale is made. While this is a reasonable assumption for “nor- 
mal" star-forming galaxies, we showed in Section [4.3] that it is not in other regimes. The 
scope of analysis possible with Mp, values predicted from SFRs is also limited, because of 
these in-built assumptions. 


4.4.2.2. Inferring Myr from scaling relations. Scaling relations have also long been 
used to estimate Mgr, with the community initially focussed on environmental studies at 
z = 0 via the HI deficiency parameter, using offsets from the scaling relations to identify 
systems affected by ram pressure or other environmental processes. As mentioned above, HI 
deficiency was originally defined in terms of optical diameters, but over time other param- 
eterizations based on colour and structural properties, broad-band luminosity, stellar mass 
or specific angular momentum have been used to various degrees of success (see 
[et a. [2021]for a critical comparison and a discussion of the limitations of these different ap- 
proaches). More recently, *photometric gas fractions" estimates have been used to increase 
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the size of z — 0 samples by taking advantage of large multiwavelength surveys in the UV 


and IR regimes (e.g., 2004| Zhang et al./2009 2015). 


In the context of this review, and to compare with the above H» results, Fig. shows 
that NUV-r color is the best single-parameter predictor of far, with a scaling of: 


log( Mg; /M.) = (—0.52 + 0.05)(NUV — r) + (0.95 + 0.02). 6. 


The scatter of H1 detections around this best-fit relation is however significant, at 0.39 
dex. We also caution that this relation is best used when NUV-r« 5, the regime where 
observations provide the best constraints. Other HI gas fraction predictors in the literature 
adopt linear combinations of properties such as a colour (NUV-r or g — r) and psx to 
correct for a second-order dependence on morphology that contributes to this dispersion 
(eg. 
interesting work, 2017) explore the application of (non-linear) artificial 
neural networks to estimate HI gas mass fractions from SDSS properties, based on a training 


, finding scatters of order of ~0.3 dex. In an 


set of galaxies with ALFALFA detections or upper limits. This technique improves on linear 
methods by yielding a lower scatter of 0.22 dex, but does not reproduce the gas fractions of 
data sets with significantly different selection such as GASS (Catinella et al.|2013) without 


additional cuts, demonstrating once again that ^what you get out is what you put in". 


5. THE COLD ISM AND GALAXY EVOLUTION 


'The general purpose of galaxy evolution studies is to understand what regulates star for- 
mation, allowing galaxies to grow differently over time, resulting in the full range of galaxy 
properties observed (e.g. masses, sizes, morphologies, metallicities and kinematics). Par- 
ticularly useful are the scaling relations between these properties, as they provide strong 
constraints on the mechanisms driving galaxy evolution. For example, it is not enough to 
be able to explain how we can form both spiral and elliptical galaxies, the same physics 
must simultaneously explain why spiral galaxies tend to be star-forming while most ellipti- 
cal galaxies are not. The properties and scaling relations used to constrain galaxy evolution 
processes compare quantities such as stellar masses, colours, morphologies, metallicities, or 
environment. Our purpose here is to illustrate how adding information about the cold ISM 
clarifies the picture. 

We focus our discussion on the distribution of galaxies in the 2D parameter space formed 
by stellar mass and star formation rate. This picture has replaced the color-magnitude 
diagram as the de facto tool to represent the evolutionary state of galaxies; it highlights the 
bimodality of the galaxy population, and allows samples of galaxies to be put in their larger 
context. Figure[7|shows the distribution of the xCOLD GASS sample in the SFR-M, plane, 
with the different panels showing how the atomic and molecular-gas mass fractions as well 
as tdep(H2) vary. These are “2D” versions of the scaling relations presented in Sec. so 
unsurprisingly similar observations are readily made: amongst the star-forming galaxies, 
fur varies most with M, while fa, and taep (H2) have a stronger dependence on SSFR. 


5.1. The star formation main sequence 


The dominant feature in the SFR-M.. plane is the galaxy main sequence (MS), the tight 
relation between M, and SFR traced by star-forming galaxies. The general observational 
consensus is that the MS is ~linear up to a critical mass of Mknee, and flattens out towards 
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Figure 7: Variations in fur, fm, and tae; (H2) across the SFR-M. plane. Small black dots 
are non-detections (ie. gas fractions are below the x COLD GASS survey limit of ~ 2%), 
and the contours show the full SDSS sample at 0.02 « z « 0.05. Note that the apparent 
hard lower limit of the passive cloud is an artefect due to the sensitivity of the observations 
and the SFR calculations, with the true SFRs of the quiescent object in fact a continuum 
to zero. The dashed and dotted lines indicate the MS and the +0.4 dex scatter around it, 


respectively. Figure reproduced from (2017). 


higher masses. The MS has been observed up to z — 4, with the MS shifting towards higher 


SFRs as ~ (1+ z)? (e.g. |Lilly et a1.[2013| Whitaker et al.]2014 2016), while 


the scatter around the MS of ~ 0.2 — 0.3 dex depends only weakly on M, and z (Speagle 
Popesso otal [019]. 

If the MS were perfectly linear and scatter-free at all masses, we would conclude that 
galaxies are self-similar, with more massive systems having more gas, and therefore being 
able to sustain higher levels of star formation activity. However, explaining the observed 
shape (i.e. the low-mass slope and the flattening above Mknee) and the amount of intrinsic 
scatter requires careful consideration of what sets the amount of gas available for star 
formation under different conditions, and the efficiency of the conversion of this gas into 
stard?] 

The exact functional form and normalisation of the z ~ 0 MS is sensitive to methodology 
and selection effects, as well as the calibrations used to calculate the M. and SFR values 
(e.g. Popesso et al.[2019). We use the MS definition of|Saintonge et al.| (2016), but slightly 
revisited by using the M. and SFR values from the GSWLC-2 (Salim et al.}2018), and the 


fitting function of[Lee et al.] (2015). The best fit function we adopt, based on SDSS galaxies 
with 0.01 « z « 0.05 is: 


P jolea) ] 70718 
log SFRus [Mo yr ] = 0.412 — log | 14 | 1010.59 | í 


5.1.1. The shape of the main sequence. A common explanation for the flattening of the z ~ 


0 MS above Mknee ~ 3 x 10+? Mo is the increased importance of stellar bulges (Erfanianfar| 
2016} |Fraser-McKelvie et al.|2021). The same phenomenon also naturally explains why 


the flattening decreases in strength as redshift increases (Lee et al.12015 2020). 


6The third key feature of the observed MS, its redshift evolution and the link with the gas 


contents of star-forming galaxies across cosmic time, is reviewed in detail in |Tacconi et al.| (2020). 
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In the local Universe and on the MS, while galaxies are always rotationally-dominated 


systems (Fraser-McKelvie et al.|2021), the frequency of stellar bulges indeed increases with 
M., as shown by the higher stellar mass surface densities (Saintonge et al./2016) and higher 


bulge-to-total mass ratios (from B/T ~ 0 at M.— 10? Mo to 0.3 at M.— 10! ^ Mo; |Cook 
et al.12020). The incidence of strong bars also increases along the main sequence (Gavazzi| 
2015). 

A simple explanation for these observations would be that while a bulge contributes to 
the stellar mass, it does not (typically) contribute to the SFR, explaining the flattening 
of the MS as bulges become more dominant. However, find that when 
building a main sequence only accounting for stellar mass in the disc component of galaxies, 
the overall shape of the main sequence does not significantly change. Another factor must 
be responsible for making even the discs of massive main sequence galaxies less star-forming. 
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Figure 8: Left: Mxy,, Mur and taep(H2) as a function of stellar mass for main sequence 
galaxies. Error bars are lo confidence intervals on the position of the mean values from 
bootstrapping, while shaded areas show the lo scatter. Best fitting relations are shown 
as dashed lines and their equations given in the caption. The dashed navy line shows the 
predicted log(Mzy,) obtained by combining Eq. |7| and the mass-dependence of tdep(H2) 
on the main sequence (i.e. the dashed orange line). Right: Mean global star-formation 
efficiency (SFR/H2), My, and Mgr as a function of the offset from the main sequence, 
normalised to the value of each quantity on the MS. Also shown are the results for the 
variation of tdep(H2) with MS offset from [Tacconi et al.| (2020) (dotted line) and 
(dashed line). The diagonal gray line shows the increase of the SFR; if either Mzy,, 
Mut or taep(H2) alone were responsible for explaining A(MS), the corresponding data would 
fall on this 1:1 line. 


We can also seek a gas-centric explanation for the shape of the MS. The flattening 
would be naturally explained if, as galaxies grow more massive along the MS, their cold gas 
reservoirs do not increase proportionally (or alternatively, their star formation efficiency 
decreases), causing a reduction of the SFR per unit stellar mass. The argument has indeed 
been made that the MS is a consequence of the more fundamental relations between My, 


and M, and the tight relation between My, and SFR (Lin et al.|2019 2020). This 
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is also illustrated in Figure Sp. For main sequence galaxies, logM7, rises e 0.79 logM, up 
to Mknee ~ 1019? Mo and flattens above it (blue line and points), while £4.,(H2) increases 
monotonically as e 0.2 logM. (orange line and points). Since any molecular gas present in 
a bulge-dominated system forms stars at a lower efficiency 
[2014], the increased fraction of bulge-dominated galaxies on the MS with mass may explain 
this global increase of tae, (H2). This may therefore be the main contribution bulges make 
to the shape of the MS. However, we note that the dependence of t4., (H2) on M. continues 
below 10'°Mo, suggesting that stellar bulges are most likely not be the only/main player 
in regulating the ability of cold gas to feed star formation. 

Combining the slopes from these fits produces a shape and slope for the main sequence 
that is in agreement with other studies at low redshifts (see e.g. the compilation of 
et al.[2014). For illustrative purposes, the reverse is shown in Figure [B] the scaling of My, 
for main sequence galaxies can be accurately predicted from log SFRms and taep(H2) (navy 
dashed line). 

Having been able to ascribe most of the shape of the MS to the molecular gas contents 
of galaxies, our attention should turn to the mechanisms that regulate this gas availability. 
First, we observe that My; increases along the main sequence as 0.37 logM, (Figure [Sp 
purple line). The shallowness of this relation and the significant scatter (c = 0.47) once 
again highlight the looser link between HI and star formation. It also suggests that the 
mechanisms responsible for fuelling galaxies with gas and setting what fraction of the total 
cold gas reservoirs satisfy the requirements for molecule and star formation are fundamen- 
tally responsible for setting the shape of the MS. 

For a standard stellar-to-halo mass ratio, the mass at the knee of the MS, ~ 3x 10'°Mo 
corresponds to a halo mass of ~ 10!?Mo. Tantalisingly, this is the threshold mass found 
in simulations, where accretion of gas transitions from cold to hot mode (e.g. 
2009). The cold mode is key to sustain the molecular gas reservoir and star formation, as 
the gas accreted through this channel is cooler, denser, and most of its mass ends up in 


the central star-forming disc (van de Voort & Schaye|2012). This picture is also consistent 


with the observation that even the disc component of MS galaxies above the knee is more 


gas-poor and less star-forming (Popesso et al.]2019 2020). 


An alternative explanation is that the internal mechanics of galaxies (also referred to as 
“secular processes” ) regulate the gas availability. For example, [Gavazzi et al. (2015) suggest 
that bars may play a crucial role in reducing star formation activity in the discs of massive 
MS galaxies. In their scenario, star-forming galaxies above the knee of the MS tend to 
experience bar instability, which funnels gas to the central region in a few dynamical times, 
triggers a nuclear starburst, the formation of a (pseudo) bulge, and subsequently a drop 
in the star formation activity. This scenario is consistent both with the higher fraction of 
bulges on the MS at high masses, and the lower gas fractions. 


5.1.2. The scatter of the main sequence. The scatter around the main sequence is a conse- 
quence of the star formation histories of individual galaxies at a given stellar mass. Under- 
standing the nature, evolution (both with mass and redshift), and the timescales associated 


with this scatter is a currently active area of research (e.g. SENTEN Speagle et al.[2014| |Matthee &| 
[Schayef2019] [Bert «t al 2021) Sherman ct a [2021 


There are (at least) two competing ideas to m the scatter about the star-formation 
main sequence. In the first, galaxies have SFHs that vary on short timescales, resulting in 
galaxies constantly moving up and down with respect to the MS. In the second, galaxies 
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have less “eventful” SFHs, resulting in their offset from the MS to persist on long timescales. 
The first model is supported by the “equilibrium” model for galaxy evolution (e.g.|Lilly et al. 
where galaxies regulate their star formation activity via gas inflows and outflows. In 
this model, a gas accretion event will be followed by an increase in star formation activity. 
The intense star formation will drive an outflow, depleting the gas reservoir further and 
forcing the SFR back down. The scatter of the MS is therefore linked to the time required 
for galaxies to return to their equilibrium state after a perturbation of their gas contents 
(Finlator & Davé]2008). Based on simulations, [Tacchella et al.|(2016] estimate that galaxies 
would oscillate around the MS on timescales of ~ 0.4tHubble (i.e. ~ 5 Gyr at z = 0, but as 
little as ~ 1 Gyr at z ~ 3). 

Other simulations however suggest that galaxies “remember” their star formation his- 
tories, with the present-day SFRs connected to the formation time and halo properties 
(Matthee & Schaye|2019). Under such a scenario, the offset of a galaxy from the MS varies 
on long timescales (~ 10 Gyr), and is related to its position in the cosmic web. Observa- 
tional evidence for this in the nearby Universe has been found by [Berti et al.] (2021), who 
report that at fixed M., galaxies above the MS are less clustered that those below. Using 
MaNGA (Bundy et al.|2015) data and reconstructing the star formation histories, Bertemes 
et al. (2021, in prep) find a similar result: the majority of galaxies at z ~ 0 are found on 
the same “side” of the MS where they have been on timescales of Gyrs, with short-term 
stochasticity only accounting for a small fraction of the scatter around the MS. 

The observation of a systematic difference in fH; and fx, between galaxies in the upper 
and lower half of the MS scatter (see Figures |7| and Bb) suggests that the scatter of the 
main sequence is gas-driven. It also argues against the rapid motion of galaxies around the 
MS, rather suggesting that MS offsets persist on timescales of ~ Gyrs. The result for Hi 
is particularly constraining; with such a long depletion time (taep(HI)~ 5 Gyr, see Fig. [c]. 
the HI reservoir is sluggish and does not respond quickly. Rapid variations in the SFHs on 
~Gyr timescales are therefore not expected to be seen in the global fg; and fy, values, 
and would rather work towards erasing the scaling between gas mass fraction and SSFR. 
Comparing all the various timescales that influence the baryon cycle and the star formation 
process, argue that since cosmic noon (z « 2) the star formation 
activity of galaxies is limited by the rate of gas accretion rather than by star formation 
efficiency. This is in line with the scenario where the offset from the main sequence is linked 
to halo properties, as those regulate gas availability and therefore star formation histories. 


5.2. Above the main sequence 


Galaxies in the local Universe can have star formation rates as high as hundreds of solar 
masses per year, putting them two orders of magnitude above the main sequence. Major 
mergers are the process most likely to significantly push galaxies well above the main se- 
quence, and indeed the highest SFRs in the local Universe are found in merging gas-rich 


galaxies at the time of coalescence of their nuclei (Sanders & Mirabel|1996 
2016). 


Major mergers at z ~ 0 are very rare events, making them minor contributors to the 
overall “star formation budget", despite the high SFRs they can generate. Far more signif- 
icant to this budget are the processes that may elevate galaxies above the main sequence 
by as little as a factor of ~ 2 — 3, but occur far more frequently, such as minor mergers, 
galaxy interactions, and bar instabilities. 
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Taking a gas-centric view, there are two possible channels to explain why dynamical 
processes such as these can elevate the SFR of galaxies: (1) they may increase the amount 
of gas available to participate in the star formation process, and (2) they may change the 
properties of the gas, resulting in an increased star formation efficiency (i.e. more stars 
formed per unit total gas mass). Figure [8p makes the case for both of these factors to be 
at play, with starbursting galaxies (defined here as those more than 0.4 dex above the MS) 
having excess H1 and Hz and shorter tgep(Hz2). 


5.2.1. Increased molecular gas mass. Dynamical processes can drag atomic gas, initially 
at large radii, into the central regions of the galaxies where it will convert to molecular 
gas under the higher ambient pressures. find significantly flatter gas- 
phase metallicity gradients in galaxies in close pairs, compared to isolated galaxies due to 
the infall of metal-poor atomic gas from the outskirts as a result of the interaction 
Rupke et al.|2010). There is also evidence for an increase in the H2/HI ratio as galaxies 
progress along a merger sequence Lisenfeld et al.12019). The cause 
of this is both an increase in the amount of molecular gas and the rapid consumption of 
the gas. For example, as systems move from the early to later stages of a merger sequence, 
their Hə mass fraction doubles while their total cold ISM mass (HI+Hz2) decreases, as a 
result of the HI begin dragged inwards, converted to Hə and then stars, or otherwise heated 


or removed from the system via outflows (Larson et al.[2016| |Georgakakis et al.12000 


5.2.2. Shorter depletion timescales. Extreme starbursts in the local Universe, observed in 
merging gas-rich systems, have H2 depletion times < 10? yr, an order of magnitude shorter 
than normal star-forming galaxies. The molecular gas in galaxies experiencing mergers is 
more centrally concentrated and pushed to higher densities (Solomon et al.[1997), conditions 
that increase the efficiency of the gas-to-stars conversion. 

In simulations, find an increase in gas mass fraction in merger- 
induced tides, which pumps compressive turbulence into the ISM (as opposed to solenoidal 
turbulence). This produces an excess of dense gas, which leads to an increase in SFR per 
unit total molecular gas mass (diffuse+dense). The simulations of [Moreno et al.| (2021) also 
show that during the interaction of two galaxies, a central increase of the dense gas fraction 
is what mostly drives the SFR enhancement. This enhancement is modest in the initial 
phase of the interaction and highest around the time of coalescence. This is supported by 
the observations of (2018), showing that star formation efficiency only increases 
significantly in galaxy-galaxy interactions at close distances (< 20 kpc) and when the two 
galaxies are similar in mass. 

Several other studies also report that taep(H2) is shorter by a factor of 2-4 in galaxies 
showing signs of M disturbances (minor mergers, galaxies with companions, galaxies 
with strong bars... (Young ct aL[IS96] Saintonge et i UI2] [Violino c al DOS). While 
this enhancement in star formation efficiency is modest, these kinds of interactions are 
far more frequent than major mergers, making them important in an integrated sense. 
However, it should be noted that not all mergers and interactions lead to enhanced star 
formation; this has been observed and simulated (e.g. [Di Matteo et al.[2007), and linked to 
the specifics of the orbits of the two interacting galaxies and their relative properties, such 
as their mass ratio and gas masses. 
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Quenching: The 
process of halting 
star formation in a 
galaxy (e.g. via 
depletion of its cold 
gas reservoir). 


Green valley galaxies: 
The population of 
galaxies with SFRs 
(or colors) 
intermediate 
between those of 
star-forming (blue) 
and quenched (red) 
galaxies. 
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5.3. Below the main sequence 


What makes galaxies “fall off’ the main sequence? All the so-called quenching mechanisms 
that have been proposed to explain the decrease of star formation in galaxies are related 
to the cold gas reservoirs, which have to be somehow depleted or prevented from forming 
new stars. Are the galaxies below the MS all gas-poor, and which cold gas phase is 
most affected? Can the decrease of SFR at fixed stellar mass be explained by reduced 
gas contents, or does it require a change in taep(H2) as well? Answering these questions 
ideally requires large, representative samples of galaxies below the MS with measurements 
of atomic and molecular gas, structural and star-formation properties. There are still 
very few studies of cold gas content of galaxies below the MS, partly due to a general 
perception that these systems are devoid of gas. While large integral-field spectroscopy 
(IFS) surveys are now underway, allowing detailed investigations of stellar populations, 
SFRs and chemical abundances on a spatially- (and spectroscopically-) resolved basis, the 
overlap between these and Hi and CO (even global) surveys is minimal, so that dedicated 
follow-ups of IFS data sets are required to measure cold gas content, especially below 
the MS. While this is a research area that will benefit tremendously from the increased 
sensitivity to small gas content afforded by the next-generation surveys with the SKA 


pathfinders and ALMA, several conclusions can be gleaned from existing work. 


5.3.1. Decreased gas content. Both far and fm, decrease steadily below the main se- 
quence. The gas content appears to decrease smoothly when galaxies move away from the 
MS at fixed stellar mass, without obvious discontinuities, until the sensitivity limits of the 
observations are reached (Figure[7), and even further below once non-detections are taken 
into consideration (Figure [8p). For example, galaxies that are between 0.4 and 1.0 dex 
below the MS have offsets of 0.41 and 0.47 dex in H1 and Hə gas fractions, respectively, 
compared to the averages on the MS, while these numbers fall to -1.10 and -1.21 dex for 
quenched systems (A(MS)« —1.0 dex). 

Although there is significant scatter in fg; and fy, below the MS, there is no evidence 
for a significant population of “passive” galaxies with HI or H» reservoirs that are compa- 
rable to those of MS galaxies. This is shown in Figure [7] as the galaxies below the main 
sequence, and especially those in the quiescent “cloud” are for the most part non-detections 
in both CO and Hi, with just a handful of points that scatter towards gas fractions that 
are comparable to those of the MS. This contradicts recent claims of the existence of a 
large population of passive, massive discs harbouring HI reservoirs as large as those of MS 
galaxies (Zhang et al. but see[Cortese et al.[2020). This is not to say that there cannot 
be a small minority of unusual systems that are outliers in some of these scaling relations. 
Good examples of these are the so-called Hi-excess galaxies, which have unexpectedly large 
Hi reservoirs for their low SSFRs [2018]. Radio interferometric and 
optical follow-up observations indicate multiple possible reasons for their HI excess, such 
as an external origin of the gas or unusually extended HI disks with large specific angular 
momentum, although not all these systems are well understood at this point 
pois) 

A class of galaxies that is typically found below the MS is that of the early-types, 
ie. ellipticals and lenticulars. Early studies reported low detection rates of cold gas in 


these galaxies (e.g., Gallagher et al. 


Driel & van Woerden|1991), reinforcing the perception that, although not always devoid 


Saintonge & Catinella 


of gas, these systems are generally gas-poor. Notably, the Infrared Astronomical Satellite 


(IRAS) detected large fractions of early-type galaxies at 60 and 100 um (Jura]1986 
1989), implying the presence of a cold ISM in typical galaxies in this class and thus 
motivating more sensitive observations of their cold gas content (e.g.,|Sage & Wrobel|1989 


[Roberts et al.|[1991). These observations are now superseded 
by modern, more representative cold gas surveys of early-type galaxies selected from IFS 
surveys such as SAURON (de Zeeuw et al.|[2002), ATLAS®? (Cappellari et al. 2011) and 
MASSIVE (Ma et al.|2014). In particular the ATLAS?P Hr survey by (2012) 


detected 21 cm emission from 32% of the 166 targeted galaxies, with the rate increasing 
to ~40% outside the Virgo cluster. Interestingly, a fraction of these systems were found 
to be as Hr-rich as spiral galaxies, but lacking the high column density gas typical of the 


latter. Detection rates of molecular gas were of order ~25% overall 
[Welch et al.[2010] [Young et al.[2011), with a dependence on kinematic bulge fraction (i.e., 
systems characterised by higher stellar velocity dispersions at fixed mass are less likely to 
have detectable molecular gas) rather than local environment (Davis et al.[2019). Naturally, 
detection rates depend on the sensitivity of the observations (as well as sample selection) 
and might be misleading. Thus, we compare ATLAS?P Hi and Hy measurements with 
xGASS and xCOLD GASS data in Figure [p] As can be seen, the populations of ATLAS?P 
early-type galaxies have significantly lower cold gas fractions compared to MS systems, with 
only a handful of exceptions, but those have SFRs putting them on the MS (see Figure Ph). 
In other words, statements that early-type galaxies are gas-poor systems and that they host 
significant cold gas reservoirs are not in contradiction, when gas poor is meant in comparison 
to galaxies with similar stellar mass on the MS and is not taken as synonym of devoid of 


gas. 
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Figure 9: Distribution in the SFR-M. plane (panel a) and comparison of fa, (panel b) and 
fur (panel c) between the ATLAS?P (red) and xGASS/xCOLD GASS (gray) datasets. In 
both samples, gas detections and upper limits are indicated by small circles and downward 
arrows, respectively. Blue lines in all panels show the best fitting relation mass for main 
sequence galaxies (|A(MS)| « 0.4). 


5.3.2. Longer depletion timescales. There is compelling evidence that changes in gas mass 
fraction alone cannot account for the full decrease of SSFR, but that it is accompanied 


by longer depletion times (e.g. 2012} |2016} |Colombo et al.|2020). Longer 


depletion times below the MS might be caused by an increased stability of the molecular 
gas against fragmentation in bulge-dominated systems. Evidence for this comes for instance 
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from the ATLAS?” early-type galaxies, which have longer taep(H2) than MS galaxies, due 
to increased shear in the inner parts of the galaxy (where the rotation curve is quickly 
rising and most H» is found; Davis et al.[2014) and/or increased velocity dispersion of the 
molecular gas (Dey et al.|2019). 'These longer depletion times might also be due to the 
fact that the CO(1-0) line traces both diffuse and dense molecular gas (see Section [5.5]. In 
other words, increasing diffuse-to-dense molecular gas mass ratios below the MS would lead 
to longer CO-based depletion times, even if the depletion time of the dense, star-forming 
gas was the same. As mentioned above, the Hr-emitting gas in these systems typically has 
lower column densities than in star-forming spirals, supporting the idea that the physical 
conditions of the cold gas reservoirs might be less conducive to star formation below the 
MS. 

Resolved studies that compare cold gas and star-formation properties of galaxies on 
and below the MS have the potential to shed light on the physical processes that drive 
the global trends discussed here. Several IFS studies have shown that the SFR density 
is suppressed at all radii in green valley galaxies, but especially in their central regions, 


suggestive of inside-out quenching (e.g., Belfiore et a1.|2018| [Medling et al.|/2018). In a 
recent work, (2020) found that spatially resolved H2 gas fraction (£ m, /3.) 


and SFE (Xsrn/Xm,) are suppressed compared to typical star-forming regions, and that 
both contribute to reducing the resolved SSFR in seven green valley galaxies. In an inter- 
esting complementary approach, (2021a) compared star-forming and retired 
spaxels of eight galaxies from the ALMaQUEST (Lin et al.|[2020) survey, and concluded 
that the retired spaxels are typically concentrated in the inner regions of galaxies and have 
significantly lower molecular gas fractions than the star-forming spaxels in the same galaxy, 
suggesting that quenching proceeds inside-out, at least for central galaxies. Although based 
on small samples, such resolved studies are in qualitative agreement with the findings based 
on global properties. The consensus is that a reduction of the molecular gas reservoir is a 
necessary but not sufficient condition to suppress star formation below the MS, as it must 
be accompanied by longer depletion times, as illustrated in Figure [sp. 


5.3.3. Multiple paths below the MS: the role of environment and AGN. The large scatter in 
both gas fractions and SFRs below the MS suggests that there is not a single physical process 
that is responsible for reducing the cold gas content and quenching the star formation in 
galaxies, and that there are likely to be multiple evolutionary paths through this region of 
parameter space. Indeed, we already know at least some of these different paths. 

For satellite galaxies environmental effects are important, and there is now evidence that 
active stripping of the cold gas reservoirs is ubiquitous in both galaxy clusters and large 
groups — although the overall picture across dark matter halo and stellar mass dimensions 
is a complex one, as recently reviewed by [Cortese et al.] (2021). Environmental effects could 
also contribute to suppress star formation in the central galaxies of groups and clusters, by 
preventing their access to gas replenishment from large-scale filaments. 

For central isolated galaxies, intrinsic or mass quenching (Peng et al. [2010), ie. halting 
of star formation caused by internal processes such as stellar or AGN feedback, is generally 
assumed to be more relevant. Super-winds and jets driven by starburst or AGN activity 
could impact star formation by directly removing part of the cold gas reservoir, or by 
injecting energy and turbulence into the gas and preventing cooling (negative feedback); 
locally, the same processes could also lead to enhanced star formation by compressing the 
gas (positive feedback). Numerical simulations and semi-analytic models of galaxy evolution 
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have been invoking AGN feedback as a key process to quench the star formation in the most 
massive galaxies for over two decades (e.g., 
[2005]. Indirect observational evidence for the impact of AGN activity on the ISM of galaxies 
comes for instance from the presence of large X-ray cavities in the hot ISM of the central 


galaxies of so-called cool core clusters (see e.g. review by 2012), as well as from 
the two large gamma-ray bubbles extending 50 degrees above and below the center of our 


own Milky Way (Su et al./2010 Guo & Mathews|2012). Excitingly, fast, massive atomic 


and molecular gas outflows have also been observed in massive galaxies hosting AGN (e.g. 
the global cold gas reservoirs of galaxies harbouring an AGN do not seem to be significantly 
affected, as several studies reported similar molecular gas fractions and depletion times for 
systems with and without AGN located in the same region of the SFR-M. plane 


2017 2018 2021), a result that holds even for powerful quasars 


2020). The same result was found when comparing HI gas mass fractions (e.g., 


Fabello et al.|2011b| |Ellison et al.|2019). Evidence therefore suggests that 


the biggest role played by AGN in the nearby Universe, at least when considering the large- 
scale properties of galaxies, is to keep halos hot. This throttles the accretion of new gas, 


pushing massive galaxies off the main sequence and preventing already quenched galaxies 
from resuming star formation activity (Best et al.|2005), but work continues to determine 
the impact of AGN activity on the cold ISM, especially on small scales. 


5.4. The mass-metallicity relation 


The scaling relation between stellar mass and gas-phase metallicity (i.e. the mass-metallicity 
relation, or MZR) provides additional clues to the importance of gas, and in particular the 
cold ISM, in regulating the growth of galaxies. The relation has been observed locally 


(e.g. 1979 2004) but also up to z > 2 (eg. 
2006), and down in the regime of dwarf galaxies (e.g. 2006). While stellar 


metallicities measure the integrated chemical evolution of galaxies over their entire star 
formation histories, gas-phase metallicity is a much more "instantaneous" quantity, relying 
on the interplay of accretion of (mostly metal-poor) gas, star formation, and outflows of 
(mostly metal-enriched) material. Given its deep connection with all aspects of the baryon 
cycle, the MZR has huge value in constraining galaxy formation models, both analytic and 
numerical (e& [Dave et al [2012 ii et al [9013] Ma et aO). 

The MZR has been extended to add a third parameter, in an effort to explain the 
correlated scatter of the relation. For example at fixed M., galaxies with higher SFRs are 


found to have lower gas-phase metallicities (Ellison et al.|2008). The SFR has been added 


to the MZR to form the M,-Z-SFR relation, referred to as the Fundamental Metallicity 
Relation (FMR, e.g. [Mannucci et al.|2010). There is however evidence that gas mass is 
the more fundamental parameter associated with the scatter of the MZR, with a strong 
observed anti-correlation found between My; and gas-phase metallicity at fixed stellar 
mass (Bothell al BOIS [Brown et al 2013) 
(2016) also report that the scatter of the MZR correlates with Mz,. These 
studies all find that the scatter of the MZR depends more strongly on gas mass than on 
SFR, which would make the FMR a consequence of this more fundamental relation between 
stellar mass, metallicity and gas mass (via the link between gas and star formation). 
'There are important practical implications for this result. For example, the suggestion 
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of a redshift-independent relation between metallicity and the stellar mass-to-gas mass ratio 
with small intrinsic scatter (Zahid et al.|2014) might provide an alternative route to predict- 
ing gas masses from optical spectroscopy (see Sec. 2-3). The MZR also offers opportunities 
to constrain important elements of the baryon cycle that are otherwise difficult to directly 
observe For example, strong outliers from the MZR with low metallicities given their mass, 
could represent objects having recently undergone an accretion event (e.g. |Peeples et al. 


2009). 


5.5. The multiple scales of star formation: from halos to molecular clouds 


In the sections above we have highlighted how the star formation activity and metallicity of 
galaxies are tightly linked to their cold ISM properties. The H1 and H2 contents of galaxies 
vary smoothly as a function of both M, and SFR, as does taep(H2), and as a result the 
position of galaxies in the SFR-M, plane can be explained by (1) the total mass of their cold 
ISM reservoirs (Mgas =My1+Mu,), (2) the fraction of this total gas that is in molecular 
form (Mz, /Mgas), and (3) the timescale associated with the conversion of this molecular 
gas into stars (tdep, H2): 


My2 
Mgas Moss 


SFR = 


tdep,H2 


where each term on the right-hand-side is a function of M. and A(MS). As illustrated in 

Figures[7]and [5] the relative importance of the different terms on the right hand side of Eq. 
8 varies for different galaxy sub-populations. For example, the shape and scatter of the MS 
is mostly driven by Mp3, high SFRs in starburst galaxies are the result of both an increase 
in fm, (though a more efficient atomic-to-molecular conversion) and shorter t4,, (H2), while 
low SFRs in quenched galaxies are mostly driven by low Mgas. 
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Figure 10: Illustration of the connections between the cold ISM and the other components 
of the baryon cycle, focusing on the path of gas from the large scale environment all the 
way to star formation. We highlight the range of scales characteristic of each step of the 
process, and the physical mechanisms at play. The HI-to-H2 transition mostly depends on 
midplane gas pressure (e.g. [Leroy et al.[2008); with contributions 
from the radiation field (e.g. 1993), while further collapse and fragmentation of 
GMCs is linked to shocks (from SNe, for example), turbulence, and magnetic fields. 
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Interestingly, the terms on the right hand side of Equation 8 are each regulated by phys- 
ical mechanisms operating on very different scales, placing the cold ISM at the crossroads. 
This is represented schematically in Figure[I0] For example, both observations and simula- 
tions suggest a deep connection between Mgas and the angular momentum and formation 


histories of the dark matter halos in which the galaxies reside (e.g. 2001 
Obreschkow et al.|2016| [Guo et al.|2017] [Lutz et al.|[2018] [Obuljen et al.|[2019 
[Piña et al.[2021). This is an active area of research that will benefit from the sensitivity 


of the next generation of observations, especially to push the study of the galaxy-halo con- 
nection to lower mass galaxies. Such a paradigm may enable the use of HI observations to 
(indirectly) probe the very largest scales of the galactic ecosystem, the dark matter halos, 
and their histories. 

The second term of Equation 8, the molecular ratio (Ma, /Mgas), is on the other hand 
controlled by the conditions within galactic discs, in particular the gas density, strength 


of the UV field, and metallicity (e.g. 1993| |Gnedin & Kravtsov||2011). For 


example, many observations show that the atomic-to-molecular ratio correlates with the 
midplane hydrostatic gas pressure (e.g. [Blitz & Rosolowsky|2004] {Leroy et al.[2008). A very 
significant development over the past decade has been the realisation that the properties 
of GMCs themselves also depend on the large-scale properties of the galactic disc in which 


), responding to the pressure of 
2018). 
the previous picture of “universal” GMCs with a constant depletion time (e.g. |Leroy et al. 
[2008] [Bolatto et al.[2008], spurred by the observation of the systematic trends in taep (H2) 
across the local galaxy population (see Section|4.3), and the ability to map CO at resolution 
of tens of parsec in galaxies spanning a wider range of properties (e.g. |Schinnerer et al. 
Schruba et al.|2019). 

Once it is recognised that the properties of GMCs can vary within galaxies, and from 


they are found (e.g. 


their environment (e.g. This is a significant shift from 


galaxy-to-galaxy, the observed variations in taep(H2) can be partly explained by the specifics 
of the molecular gas tracer employed, the CO(1-0) emission line. With a critical density 
of ~ 10? cm ?, a significant (but variable) fraction of the total CO(1-0) emission traces 
diffuse gas, while star formation mostly takes place in the densest regions of molecular 


clouds (Heiderman et al.|2010). In the Milky Way, for example, 2596 of the molecular gas 


is diffuse (defined as gas detected in !?CO but not in !* CO), in fact only about 15% of the 
molecular gas mass traced by !?CO is in observed molecular cloud complexes 
et al.]2016). 

Variations due to changes in the dense-to-diffuse molecular gas ratio can be captured by 
explicitly distinguishing between the total My, (as traced by e.g. CO(1-0)) and the dense 
star-forming gas in Equation [5.5} 


. Maz , Maense 
gas ` Moo. ^^ Mao 


SFR = 


, 9. 
tdep,dense 

where Maense represents the mass of dense molecular gas, say with density in excess 

of ~ 10° cm °, and taep,dense = Maense/ SFR. Efforts are ongoing to determine how 

Maense/ M2 and tdep,dense might vary systematically both within galaxies and from galaxy- 

to-galaxy, so far with a focus on using HCN and HCO*. With higher effective critical 

densities even for their low-J lines, these molecules probe the star-forming gas more directly. 


Gao & Solomon) (2004) indeed reported a linear relation between Lyon and Lir (the latter 


a proxy for SFR) and argue that SFR depends on the amount of dense molecular gas traced 
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by HCN emission, rather than the total molecular gas mass traced by CO emission. 

The EMPIRE program was designed to look for systematic varia- 
tions in the dense gas ratio and tdep,dense by mapping HCN(1-0), HCO? (1-0), HNC(1-0) 
and CO isotopologues at ~1-2 kpc resolution in nine nearby massive spiral galaxies. The 
observations reveal that the HCN/CO ratio varies systematically within galactic discs and 
is highest in high pressure environments (Usero et al.|2015), as shown in Figure (left). 
Interestingly, even after accounting for the different conditions within the discs of each 
galaxy, galaxy-to-galaxy variations in the HCN/CO ratio remain, correlating with global 
galaxy properties such as jj. and the Hr-to-H» ratio (Jiménez-Donaire et al.[2019). 

There are also systematic variations in star formation efficiency with global galaxy 
properties, even when only accounting for the HCN-traced dense gas (Figure[11] right). This 
result is particularly interesting, as it goes against simple expectations, even showing that 
the dense gas star formation efficiency is lower in high pressure regions. This could be telling 
us that HCN does not perfectly trace the dense star-forming gas, and that with a different 
tracer with an even higher critical density we would reach the point of a constant tdep,dense 
as proposed in certain theories of star formation (e.g. [Evans et al.|2014]. Alternatively, 
Figure[11] (right) could instead be telling us that our expectation of a universal dense gas 
depletion time above a certain gas density threshold is not correct, which would argue 


in favour of turbulence-driven theories of star formation (e.g. |Federrath & Klessen||2012). 


This is challenging work with current facilities. For example, HCN(1-0) is ~15-30 times 


fainter than CO(1-0) emission (Gao & Solomon|2004), limiting studies to modest samples 


of mostly main-sequence and starburst galaxies. Further studies will however be key to 
quantify scaling relations and test star formation theories, as discussed in Section [6] below. 
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Figure 11: Variations in the dense-to-diffuse molecular gas ratio as measured by the 
HCN/CO luminosity ratio (left) and in the star formation rate per unit dense gas mass 
(right) as a function of local pressure in 9 galaxies from the EMPIRE survey. Figure repro- 


duced with permission from |Jiménez-Donaire et al.| (2019). 


6. CLOSING THOUGHTS AND OUTLOOK 


'The past couple of decades have seen tremendous progress in the study of the cold gas in 
galaxies, with the undertaking of several large, dedicated observing campaigns that have 
measured the atomic and molecular hydrogen content of sizeable samples of galaxies on 
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global and local scales, and placed these in the context of their optical and star-formation 
properties, as well as environments. In this review we have focused on global H1 and H2 
scaling relations in the local Universe (z < 0.05), and discussed the link between cold gas 
content and star formation of galaxies across the SFR- M. plane. 

Overall, the empirical global cold gas scaling relations seem well established, at least 
for galaxies with stellar masses M. > 10°Mo, for which the systematic variation of atomic 
and molecular gas fractions has been quantified as a function of stellar mass and other 
properties (e.g. Figure [4). For galaxies with stellar masses below this threshold, global Hr 


measurements do exist (e.g., 2015), but systematic CO 
measurements are still rare and dominated by upper limits (e.g., 
eet al.]2017). These have been compiled and homogenised by (2018) | 
to extend the cold gas fraction-stellar mass relations down to M. ~ 10° Mo. 
There is tantalizing evidence for a break at M. ~ 10°Mo in the Hi scaling relation (see also 


Maddox et al.|2015| based on ALFALFA), but more statistics is needed to confirm these 
results. 


Furthermore, much remains to be learned about the physical drivers of these relations 
and how the global gas properties of galaxies are linked to their molecular clouds and 
dense core distributions, and on the other side to the large scales via the surrounding CGM 
(which extends roughly from the stellar disk to the virial radius of a galaxy), where inflows 
and outflows contribute to replenish or reducing the cold gas in the ISM. In what follows, 
we briefly reflect on how our view of galaxies and their gas reservoirs has changed in the 
past decade and highlight a few areas where future work will allow us to make progress in 
this field. 


A changing view of the gas reservoirs of galaxies Galaxies are no longer viewed as 
isolated entities, but as systems embedded in a cosmic web whose evolution is regulated by 
exchanges of gas with the surrounding environment. This gas flows from the Intergalactic 
Medium (IGM) through the CGM and into the ISM, where it is processed into stars, and 
then partly returned to the CGM/IGM via outflows from supernova explosions, stellar 
winds, or nuclear activity. Through this cycle, properties such as temperature and density 
of the gas change by orders of magnitude, from the tenuous, hot ionised IGM to the dense, 
cold neutral cores of GMCs. Focussing on the ISM (the topic of this review), various 
tracers are used to measure the global gas content in its different physical states, each 
with its associated depletion time, which is the time needed to consume the total amount 
of gas in the given phase via star formation, assuming constant SFR and no recycling. 
The fact that depletion times increase when moving from dense to diffuse gas phases 
means that they somehow encapsulate information on how close the given gas phase is to 
the star-formation process — in a way, akin to clocks that count down to star formation. 
Although the inverse of a depletion time is usually referred to as star-formation efficiency, 
this quantity tells us very little about how star formation proceeds at the scale of GMCs, 
where new stars are actually formed (efficiency should also be a dimensionless quantity). 
In this context, all these gas phases are transient, i.e. way stations toward star formation, 
rather than “reservoirs” — although it is convenient to picture these as such when depletion 
times are of order of a Gyr or more. In order for a galaxy to form stars, its dense gas phase 
needs to be constantly re-supplied, ultimately from the ionised IGM, which is the reservoir 
that fuels star formation at all epochs and is replenished by outflows. 
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Connecting ISM and CGM. The interface between the ISM and the CGM, which is 
where gas is exchanged and recycled, is of crucial importance for our understanding of 
galaxy evolution, but we lack an empirically motivated model for how these gas flows 
proceed (see however [Gaspari et al.|[2013], and how they may vary across the galaxy 
population. This is an important topic beyond the scope of this review, but it is nonetheless 
worth identifying areas where future cold gas observations might help. Because of its low 
density the CGM is best studied via quasar absorption spectroscopy, a field that has 
benefited enormously from the availability of the Cosmic Origins Spectrograph on the 
Hubble Space Telescope (Tumlinson et al.|2017). Interestingly, studies of galaxy-quasar 
pairs (where galaxies happen to be near the line-of-sights to background quasars) are 
helping to shed light on the connection between ISM and CGM, by investigating how 
the strengths of ultraviolet absorption features tracing the CGM vary with impact 
parameter and galaxy properties such as HI content and SFR (e.g. [Borthakur et al] 


2013) |2015; |2016). Complementary to this approach, surveys probing multiple quasar 


sightlines around the same nearby galaxy, such as project AMIGA (Absorption Maps 
In the Gas of Andromeda; [Lehner et al.][2020), can better illustrate how the structure 
and properties of the CGM vary with galactocentric distance, perhaps in relation to 
a more direct impact of galactic feedback in the inner regions of the CGM. These 
observational efforts are accompanied by significant efforts to simulate the CGM and 
its connection with the ISM at high resolution (e.g. 
fet _al.][2019). In the future, larger samples of galaxies with resolved and sensitive HI 
and H5 measurements, probing the low column density outer regions of galaxy disks, 
will allow us to extend these studies and explore the ISM/CGM connection in greater detail. 


Connecting global gas and local star formation. Further down the line of the baryon 
cycle, tremendous progress has also been made over the past decade to connect the process 
of star formation to global galaxy properties. There is nonetheless still a striking disconnect 
between our ability to investigate the link between gas and star formation at the scale of 
individual cores in the Milky Way, and studies of external galaxies, where we are not yet 
able to probe the spatial scales at which the physics of star formation is important. A key 
challenge for galaxy evolution studies is to connect these widely different scales across the 
galaxy population. Surveys of cold gas in nearby galaxies at kpc resolution (e.g., THINGS, 
HERACLES) have been invaluable to start bridging this gap, and high-resolution studies of 
individual galaxies at tens of pc resolution (e.g. PAWS, |Schinnerer et al.]2013) are pushing 
this even further. The PHANGS survey has carried out arcsecond CO(2- 
1) imaging of 70 nearby star-forming galaxies with ALMA, reaching the size of typical GMCs 
(7100 pc resolution). These efforts are providing key empirical constraints on the physical 
link between star formation and gas near the “cloud” scale and the galaxy-scale environment 
(o. 
2021). Further progress will require even larger samples, extending the 
parameter space especially towards low mass galaxies and below the main sequence. The 


synergy between high-resolution, sensitive cold gas studies and IFS surveys will be critical 
to advance this field. 

'The connection between gas and star formation will further benefit from exploiting a 
larger range of cold ISM tracers, in particular molecular lines with high critical densities 
such as HCN, HCO*, CS, or even higher-J CO lines. To expand on the current state of the 


art (e.g. the EMPIRE (Jiménez-Donaire et al.]2019) and MALATANG (Tan et al.|/2018) 
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surveys), we will need larger samples, especially including galaxies below the star-forming 
sequence where CO-based depletion times are longer. Even if spatial resolution has to be 
sacrificed to assemble such sample (or if spectral stacking has to be used), there would 
be much to be gained from dense gas versions of the scaling relations presented in this 
review. Comparing the shape and scatter of the scaling relations from diffuse and dense 
gas tracers will shed light on the role of secular and environmental processes regulating the 
global star formation activity of galaxies. Mapping the dense gas at resolutions of tens of 
pc in galaxies displaying a wide range of internal conditions will also be crucial to put to 
the test various models for star formation, for example confronting density threshold and 
turbulence-driven models. 


In conclusion, much progress is being made towards our understanding of the gas- 
star formation cycle via multi-wavelength, multi-scale investigations. For global gas mass 
measurements, further breakthroughs will require large representative samples that ex- 
tend the current parameter space coverage, especially pushing into the regime of low mass 
(dwarf) galaxies (M. « 10? M). The Five-hundred-meter Aperture Spherical radio Tele- 


scope (FAST; will have a role to play in these efforts for HI studies, 
complementary to that of wide-area, blind Ht surveys such as WALLABY (Koribalski et al. 
2020) on the Australian Square Kilometre Array Pathfinder (ASKAP [oT 
eae Similarly, the next generation of (sub)millimeter single-dish telescopes 
(for example the AtLAST concept?) will allow for orders of magnitude increase on the cur- 
rent molecular gas sample sizes. However, extending large molecular gas surveys into the 
regime of low mass (and therefore low metallicity) galaxies will first require further cali- 
bration of our tracers of the molecular ISM, as both CO- and dust-based methods present 
challenges under such conditions. A next generation far-infrared telescope, such as the 
Origins Space Telescopd?] would significantly help in these efforts. 

Our view of the link between gas, star formation, and global galaxy properties will be 
revolutionised by resolving atomic and molecular gas in large samples containing not only 
star-forming galaxies, but spanning the full range of stellar masses, environments, size 
and angular momentum, and including systems lying significantly off the main sequence 
(either above or below). This has been particularly challenging for cold gas surveys so far, 
traditionally limited by a trade-off between spatial resolution and number statistics, and 
for Hr emission studies also restricted to the very local Universe for individual detections 
(i.e, z < 0.4; [Fernández et al.|[2016). But there is much to 
look forward to in the coming decade(s), thanks to the powerful synergy between IFS 
instruments on 10 meter-class (and higher) telescopes on one hand, and ALMA, the Square 


Kilometre Array (SKA; 2009) and its precursor facilities (e.g., ASKAP 


and the South African Meer-Karoo Array Telescope, MeerKAT, |Jonas & MeerK AT 
2016) on the other. While the former will measure star-formation, metallicity and 


stellar properties in exquisite detail, the latter will probe the cold gas component with 
unprecedented sensitivity and resolution, also pushing Hı detections to redshift 1 and 
beyond. 


Thttps://www.atlast.uio.no 
Shttps://asd.gsfc.nasa.gov/firs/ 
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SUMMARY POINTS 


iL. 


The cold ISM in the nearby Universe accounts for ~ 1.3% of all the baryons, and 
35% of the total mass density of stars. Massive galaxies (M. 10° Mo) account of 
~ 65% of all this molecular gas, but only ~ 40% of the atomic gas. 

The atomic and molecular gas mass fractions vary systematically across the galaxy 
population. There is a sharp decline in both atomic and molecular gas contents 
above M,> 10'°?Mo. Below this threshold, the My,-to-My ratio decreases 
steadily; low mass galaxies have extended atomic gas reservoirs that are largely 
decoupled from the star formation process. 

The depletion timescale for HI is on average ~3 Gyr, but with very large galaxy-to- 
galaxy variations which dominate over any systematic trends. On the other hand, 
tdep(H2) varies systematically, in particular as a function of SSFR (or offset from 
the main sequence). These variations are likely caused by changes in the dense-to- 
diffuse molecular gas ratio. 

The shape of the star formation main sequence is set mostly by the dependence of 
My, on M. and to second order by a slowly increasing taep(H2) with mass. Star- 
forming galaxies just above (below) the main sequence have higher (lower) molecular 
gas fractions, suggesting that the scatter of the main sequence is gas-driven and that 
main sequence offsets persist on timescales of ~Gyr. 

Starburst galaxies reach their high levels of star formation due to a combination 
of increased molecular gas content (via a more efficient atomic-to-molecular gas 
conversion) and decreased global depletion times. 

Below the MS, galaxies are mostly characterised by lower cold gas reservoirs, with 
a second-order contribution from longer depletion times (again, likely related to 
variations in the dense-to-diffuse molecular gas ratio). There is no evidence for a 
significant population of quiescent but gas-rich galaxies. 

Overall, the position of galaxies in the SFR-M.. plane can be explained by their 
gas contents and the global efficiency of the gas-to-star conversion process. This 
highlights the importance of phenomena on a huge range of physical scales in driving 
galactic evolution: from cosmic web-regulated gas accretion, to cloud-scale star 
formation. 


FUTURE ISSUES 


1. 


Extending current efforts to establish cold ISM scaling relations into the regime of 
low mass galaxies (M.< 10? Mc) beyond current studies is needed to better under- 
stand the baryon cycle and processes such as cosmic reionisation and the galaxy-halo 
connection. This will require the sensitivity of new facilities, and significant work 
in understanding the systematics of the cold ISM tracers in extreme environments 
(e.g. very low metallicities, or very strong UV radiation fields). 

Mapping of atomic and molecular gas (using both diffuse and dense gas tracers) at 
high resolution (tens of pc) across a broad range of galaxy properties and well into 
the quiescent population will be required to uncover the multi-scale nature of the 
star formation process. 
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3. Exploring the connection between the cold ISM and the other gaseous phases of 
the galactic environment (in particular the CGM and IGM) as well as their dark 
matter halos will be key in understanding what regulates the availability of gas and 
therefore, down the line and many orders of magnitudes in physical scales, the star 
formation activity of galaxies. 
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